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We have sampled the diffuse 157*74 pm [CII] emission from the 
Galaxy at several galactic longitudes near the galactic plane including 
complete scans across the plane at A 11 ** 2.16° and A* 1 = 7.28°. The 
observed [CII] emission profiles follow closely the nearby 12 C0 (J»l+0) 
emission profiles. The [CII] emission probably arises in neutral photo- 
dissociation regions near the edges of giant molecular clouds (GMC's). 
These regions have densities - 350 cm*" 3 and temperatures T ~ 300°K, 
and amount to w 4*10 8 Mq of hydrogen in the inner Galaxy* We estimate 
the total 157.74 pm luminosity of the Galaxy to be ~ 6*10 7 Lg, The 
volume filling factor of the [CII] emitting regions is ~ 10~ 3 . 

We have also made estimates of the galactic emission in other far- 
infrared (FIR) cooling lines. We find the [CII] line to be the dominant 
FIR emission line from the galaxy and the primary coolant for the warm 
neutral gas near the galactic plane. Other FIR cooling lines predicted 
to be prominent in the galactic spectrum include the 51.82 pm and 88.36 
pm fine structure lines of 01 II and the 63.17 pm and 145,53 pm [01] fine 
structure lines. The OIII lines are principal coolants for the ionized 
gas while the 01 lines are significant coolants for the warm neutral 


We have also measured the 145,53 pm [01] emission line from the 
Orion Nebula with both a rxX* beam centered on the Trapezium and a 

a a 

7'x7' beam encompassing most of the Orion Nebula* By comparing our 
measured 145,53 pm intensity with previous measurements of the 63,17 pm 
COI] and 157,74 pm [CII] intensities we determine that these three cool- 
ing lines all arise in the warm (T ~ 350°K) dense (n^ ~ 2,3x10 s 
cm" 3 ) carbon photodissociation region just beyond the ionization front 
separating the compact; HI! region from the molecular cloud. We find the 
63,17 pm radiation to be self-absorbed with line center optical depth 
t ~ 2 and Doppler width (full -width-half-maximum) ~ 3,2 km s m1 . 

We have made these measurements using both NASA's Lear Jet facility 
and Kuiper Airborne Observatory (KAO), The Lear observations were made 
using a liquid helium cooled grating spectrometer which has been fitted 
with stressed Ge:Ga photoconductive detectors. This instrument is sen- 
sitive in the spectral range 100 pm < \ < 195 pm with a spectral resolu- 
tion R s \/&\ - 200, For the KAO observations, we use a modified 
version of our grating spectrometer as an order sorter for our lamellar 

grating interferometer, with which we have achieved spectral resolutions 

* 

in excess of 2000. 

The analysis of the Cell] and [01] lines presented in this thesis 
is a more sophisticated treatment of previous work which has appeared in 
the Astrophysics! Journal (Stacey et al. 1983a, b and 1985). 
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"We changed again, and yet again, and now it was too late 
and too far to go back, and I went on. And the mists had 
solemnly risen now, and the whole world lay spread before me. 


Pip- in Charles Dickens' Great Expectations 


CHAPTER 1 


INTRODUCTION 


The far-inf pared (FIR) may be loosely defined as those wavelengths 
In the electromagnetic spectrum lying between 20 and 300 pm, Spectro- 
scopic work done at wavelengths beyond 300 pm, often termed "submilli- 
meter", is distinguished from FIR spectroscopy,, in that submillimeter 
radiation is usually detected with coherent techniques, analogous to 
heterodyne receivers in use in modern radio astronomy. In contrast, 
until very recently, most spectroscopy done in the FIR has been done with 
incoherent devices in combination with optical filtering techniques. 
Instruments in current use for FIR spectroscopy include grating spectro- 
meters, Fabry-Perot interferometers, Michelson interferometers and the 
lamellar grating interferometer discussed in this thesis. The various 
devices have their advantages and disadvantages which we discuss in 
Chapter 2, 

There is a wealth of astronomical information to be gained through 
studies in the FIR. The FIR is unique in several respects: 

1) Objects as different as planets and quasars emit the bulk of 
their radiant energy in the FIR, and show the peak of their spectral 
emission at wavelengths longward of 20 pm. 


2) The FIR is unhindered by the extinction problems which plague the 
optical and near Infrared astronomers. Typical grain sizes in the 
general interstellar medium (ISM) are « .1 pm so the Galaxy is 
virtually transparent to FIR radiation. Regions totally obscured to 
optical observers, such as the galactic center, are accessible to FIR 
experimenters, 

3) The FIR is filled with emission lines produced through collision” 
al excitation of atoms, ions and molecules with low-lying energy states. 
These lines are usually optically thin, and have excitation temperatures 
30° < T < 600°K. Thus, these lines are sensitive temperature probes in 
the intermediate temperature regimes associated with low-velocity shock 
or ionization fronts or warm neutral gas. The optical lines are not 
excited in such low temperature regimes, while the radio emitting levels 
are largely populated according to their statistical weights and are 
therefore insensitive to temperature, 

4) In high temperature ionized gas regions, the FIR emitting level 
populations are temperature insensitive and determined primarily by the 
densities in the emitting regions (Chapter 3), The FIR lines provide 
excellent density probes for these regions. 

Unfortunately, the Earth's atmosphere is largely opaque in the FIR 
due to the omnipresent molecule, H 2 O. However, the column density of 
water vapor declines dramatically with increasing altitude, decreasing by 
a factor of r 500 as- one ascends from mountaintop altitudes (4,2 km) to 
aircraft altitudes (14 km) (Traub and Stier 1976), Thus, one is able to 
pursue sensitive observations in the FIR from NASA's two airborne obser- 
vatories, the Lear Jet facility and Kuipef Airborne Observatory (KAO), 


Trie Lear Jet Is the smaller of these two observatories, containing a 
30 cm f/6,5 telescope mounted In an open port In the airplane's fuse* 
lage. The power Incident on detectors from an Infrared source whose 
angular extent, Q$» is greater than the telescope beam, flp is 
proportional to Aa-j. where A is the surface area of the telescope 
primary mirror. The second law of thermodynamics requires that this 
product be a constant for any diffraction limited telescope regardless of 
aperture size. Thus, for diffuse emission regions, the light gathering 
power of the Lear Jet telescope is the same as that for the much larger 
telescope aboard the KAO, Furthermore, the larger chopper throw of the 
Lear Jet permits observations of more extended regions {see Section 
2.1). Thus, for observations of extended sources of emission (e.g. 
galactic plane observations — see Chapter 4) we use the Lear Jet 
facility. 

For sources of angular extent, smaller than the telescope 
* 

beam size, Qp the power incident on our detectors is proportional to 
An<j, Thus, for sources which do not fill the 4'x7' beam of the Lear 
Jet, the 91 cm telescope aboard the KAO has up to ten times the light 
gathering power of the Lear telescope. The KAO is a much larger facility 
(being a modified C141A Starlifter) and thus may accommodate a host of 
experimenters and crew in addition to larger instrument packages on the 
telescope. Therefore, for high spatial resolution (~ l'xl 1 ) measure- 
ments, or for measurements of point sources, we use the KAO facility. 

Observations reported in this work come from both the Lear Jet and 
KAO facilities. Lear Jet observations were made with a liquid helium 
cooled grating spectrometer, which is essentially an improved version of 
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the spectrometer first described by Houck and Ward (1979). We have 
Installed stressed GesGa detectors In this system which gives us sensi- 
tivity In the spectral range 100 pm £ \ £ 195 pm. We have also improved 
system sensitivity through the implementation of new detector designs and 
manufacturing techniques and the Installation of low noise J-FETs in the 
first stage electronics. 

The additional space available aboard the KAO allows us to employ a 
lamellar grating Interferometer for our KAO observations. We use the 
grating spectrometer as an order sorter for the resolution achieving 
Interferometric stage. The spectrometer used for KAO observations Is 
Identical to that used for Lear observations except that the f/6.5 Lear 
optics have been replaced with f/17 optics to match the beam of the KAO 
(Mel nick 1981). Our instruments are described in detail in Chapter 2, 

The first astronomical detection of a FIR fine structure line was 
the observation of the 88 pm [OIII] line emanating from the M17 (Ward et 
; al. 1975). Since then, a host of FIR fine-structure and molecular 
rotational transitions have been observed in emission and absorption from 
interstellar clouds. A compendium of FIR fine structure lines is 
presented in the Appendix. 

The FIR fine structure lines arise in atoms with 1, 2, 4, or 5 equi- 
valent p electrons. The spin-orbit (L-S) interaction splits terms of 
total spin, S, and total orbital angular momentum, L, into a set of fine 
structure levels, each of which represents a different value for the 
total angular momentum vector determined by the vector addition of L 
and S. Since transitions among states of different J involve no change 
in the electronic spin or orbital angular momentum, they are forbidden to 
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electric dipole radiation. Magnetic dipole transitions, however, are 
allowed, and as such the excited states are metastable with lifetimes of 
the order of days. The excited levels may thus be colllslonally depopu- 
lated In dense (nm £ 10 s cm" 3 ) regions. We discuss the physics of 
fine structure line emission and radiative transport In Interstellar 
clouds In Chapter 3. 

Since the FIR fine-structure lines are generally optically thin and 
colllslonally excited at low densities, they serve as excellent coolants 
for the general ISM (see Section 3,3), The dominant coolants for warm (T 
> 70°K) low density (nn < 10 6 cm" 3 ) neutral clouds are the fine struc- 
ture lines of the astrophyslcally abundant species CII and 01. As col- 
lapsing clouds must go through these temperature and density ranges when 
forming protostars, studies of the [Cl 13 and [013 lines In star formation 
regions should lead to a better understanding of how stars are formed 
from interstellar gas. We have measured the 3 Pq - 3 Pi (145.53 pm) [013 
transition in the Orion Nebula, believed to be a region of active star 
formation. In Chapter 5, we Incorporate this measurement with previous 
measurements of the 3 P 3/2 - 3 Pi /2 (157.74 pm) [CII] fine structure 
line to build a model for the warm neutral regions just outside the 
Ionization front in Orion, 

The primary thrust of this work, the measurement of the 157 pm [CII] 
emission from the Galaxy is contained in Chapter 4, Most of the mass in 
the galactic interstellar medium is contained within atomic or molecular 
hydrogen clouds. Neutral clouds near the galactic plane are probably 
enveloped in warm (T - 3Q0°K) low density (nfu < 10 3 cm" 3 ) "CII" 
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regions generated by the Ionising effects of the Interstellar radiation 
field (see Section 4.3). The dominant cooling transition for these warm 
neutral regions is the [0113 fine structure line at 157.74 pm, with a 
lesser contribution from the 63.17 m E0I3 line. We have sampled diffuse 
ECU] emission from the Salary at several galactic longitudes Including 
complete scans at *** * 246° and « 7,28*. These measurements 
allow us to deduce that the ECU] emission arises predominantly from the 
edges of giant molecular clouds within the solar circle and that there 
are knotlike features In the ECU] emission. We are able to estimate 
temperatures and densities In the emitting regions and estimate the ECU] 
luminosity and the total ECU] emitting mass of the Galaxy. We also 
predict the galactic luminosity of other FIR and submllllmeter emission 
lines. We conclude that the ECU] line is the dominant FIR emission line 
In the Gal aw and the primary coolant for the warm neutral gas. 


CHAPTER 2 
INSTRUMENTS 

" * 

2.1 » .Introduction^ 

The primary design consideration for any FIR spectrometer Is the 
detection of the maximum number of astrophysically Interesting spectra] 
lines In the minimum time. Detection of these lines is hampered by the 
fact that the atmosphere, being a good absorber of Infrared photons, is 
also a copious emitter of Infrared photons. In addition to this "sky" 
background the telescope and telescope cavities are also good Infrared 
emitters, adding to the background flux, bet us Investigate the nature 
of this background. 

If v/e assume the background emitters behave like a graybody of 
emisslvity e v at temperature, T, then the total number of background 
photons at frequencies, v, absorbed by the detector per unit time, in 
frequency interval dv 1st 

N v “ tvTt v B v e vAD dv/hv 

where 

B v » Planck blackbody functions 

nv * quantum efficiency of the detector, i.e., 

# photons absorbed * # photons incident on the 
. : detector; ■■ 
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t v ■ transmission of helium cooled optics. 

An « aperture - solid angle product of the optical 
system On the Ideal, diffraction limited case 
An * 3.7 \ z cm 2 sr )| 
h * Planck's constant. 

In terms of the system resolution, R s v/dv we have 

Ny ** nytyeyByAa/hR 

For the present calculation, we assume e v Is given roughly by the 
emlsslvity of the sky at 150 m and at aircraft altitude of 41,000'; 
e v « ,25 (Traub and Stler 1976), Taking a "typical" transmission for 
a grating spectrometer t v w * 3 * detector quantum efficiency n v a 
,15, and a sky temperature T - 25U°K we have; 

N v - (4*10 n )/R photons/sec 

If we let R * 200 we will have My ~ 2*10 9 photons/sec. In con« 
trast, from an especially strong continuum source at 150 m we might 
detect ** 1G 3 photons/sec at this resolution, and from a strong line at 
this wavelength we might detect 2xl0 7 photons/sec. Thus, we are faced 
with the problem of detecting source radiation which is at least a 
factor of 50 weaker than the background radiation. The procedures 
generally employed to compensate for this large thermal background are 
called "chopping" and "nodding". The Cassegrain telescope secondary 
mirrors on both the Lear Get and KAO are designed to wobble in- a square 


wave pattern thus causing the beam alternately to look at two portions 
of the sky. If one places the FIR source In one of the two beams, the 
result will be an oscillating photon flux at the detector due to the 
source, on top of a large D.C. offset due to the thermal background. 

This procedure Is called “chopping" and the two beams are often referred 
to as "plus" beam and "minus" beam. Unfortunately, It Is nearly always 
the case that one side of the telescope cavity Is warmer than the 
other. Thus, the thermal background may be larger In one beam than 
In the other. There might also be a gradient In the sky background.' 
These gradients in the thermal background will then give rise to an 
A.C. signal. To compensate for this second order effect, one uses the 
procedure called "nodding", i.e., first Integrating with the source In 
the plus beam for a certain integration time (usually " b seconds) and 
then Integrating again with the source In the minus beam. The differ- 
ence between these integrations will then give twice the signal due to 
the source and effectively eliminate first order thermal gradients. 

The above procedures enable us to distinguish a small source signal 
in the presence of a large thermal background. Unfortunately, these 
procedures do not eliminate the noise associated with this background 
flux. 

In the most optimistic case, the noise in the background flux is r 
due to the statistical fluctuations in photon arrival rates. If we let 
N y * number of background photons absorbed by our infrared detector as 
above, then the RMS fluctuations in the photon flux is given by (Boyd 
1982) 
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ANy * C 2Ny (l+flyny)/At 3 l /2 

where 2 Is the photon occupancy and at Is the Inte- 

gration time* The first term in parentheses Is the classical "shot 
noise" due to the discrete arrival rate of photons. The second purely 
quantum mechanical term Is often referred to as "excess noise" and Is 
due to photon bunching* Due to the Bose*E1nste1n nature of the radia- 
tion field the mean number of photons per field made Is greater than 
unity. For T - 250K at 150 urn* n - 2.1. 

In terms of the system parameters resolution, R, and throughput 
T v 5 T v n v the fluctuations seen by the system at 150 urn are 

aN v - 5xlQ 6 (Ty/R ) l/2 

Thus, one may reduce the background fluctuations by Increasing 
system resolution, or decreasing system throughput. Obviously one does 
not wish to minimize the fluctuations by decreasing system throughput as 
the signal strength Is proportional to T v and the photon noise Is only 
proportional to Tj^ 2 . Therefore, the way to decrease system noise, when 
the system noise Is dominated by photon noise, is to Increase system 
resolution. . 

With recent developments In crystal growth and electrical contact- 
Implantation techniques, Infrared detectors are now available whose 
Intrinsic noise Is much less than the photon noise above. Johnson noise 
may be reduced through suitable selection of load resistors and ampli- 
fiers are also available whose intrinsic noise is much less than the 
photon noise. One should, then, be able to achieve photon-noise-llmlted 
performance in an aircraft environment. Generally, however, It is found 


u 

that background fluctuations are a factor of 5 or 10 larger than those 
due to photon noise. The cause of this excess noise Is probably a 
combination of changing airmasses, fluctuating sky emissivlty and 
temperature* or atmospheric turbulence Inside and outside of the tele- 
scope cavity. However* the conclusion remains: higher system resolu- 
tion should result In smaller background noise. For discussions of 
background-limited performance In the FIR, see Keyes and Qulst (1970), 
Glllett, Dereniak and Joyce (1977), and Boyd (1982). 

Increasing system resolution is beneficial In other respects also. 

1) Many FIR spectral lines are formed in regions which also have 
high continuum flux levels due to thermal emission from interstellar 
dust. Until one actually resolves a line, higher spectral resolution 
will enhance the presence of a line at the expense of the continuum. 

This is important as the system response Is a convolution of many 
functions which vary from one wavelength to another. Included are the 
spectrometer response function, filter transmission functions, atmo- 
spheric transmission and beam profile. Much of this problem Is allevia- 
ted through dividing the source spectrum by a suitable calibrator spec- 
trum which Is presumed to emit like a blackbody (ase Section 2.6). 
However, false spectral features may still persist. These may be caused 
by different solid angles being subtended by the source and calibrator 
tending to alter the instrumental profile; wavelength calibration 
errors due to pointing errors; different alrmasses through which the 
source and calibrator were observed or even spectral features In the 
calibrator spectrum Itself. These problems are minimized through the 
use of a high resolution spectrometer which need only observe the 
continuum immediately adjacent to the interstellar spectral feature. 
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A good example of calibration problems Is the 63.17 pm £01] line 
which Is located on the wing of a very deep telluric water vapor absorp- 
tion feature only - .15 pm away. If the resolution of the spectro- 
meter used to measure the 01 feature is only moderate ft/M ** 200), 
the water vapor and £01] features are not distinguished by the spectro- 
meter. Hence, it Is necessary to obtain a calibration spectrum whose 
wavelength scale Is precisely that of the source spectrum and In which 
the water vapor feature Is saturated to the same degree as the source 
spectrum. In practice, this Is not easy to do and large calibration 
errors may remain. If the saturation of the water feature Is higher In 
the calibrator, the divided spectrum might show an emission feature to 
appear at 63,17 pm even If the source had no such featurel Higher 
spectral resolution eliminates much of this problem. 

2) A second reason that high resolution Is desirable is that when 
an Interstellar line Is actually resolved, one may Infer from the line 
width bulk physical properties of the emitting region such as turbulent 
velocities, or If the broadening is thermal, kinetic temperature. 

3) Lastly, the FIR spectrum many regions is so thick with spec- 
tral lines, that one needs high resolution to separate and correctly 
identify spectral features. A good example of this is the BN-KL region 
in Orion where we have recently discovered the 163.124 pm and 163.397 pm 
emission Tines of excited OH (Vlscuso et al. 1985), These lines are 
immediately adjacent to the previously measured 162.812 pm, J * 1.6*15 
rotational transition of shocked CO. If the resolution of our instru- 
ment had been less than - 500, it would have been Impossible to dis- 
tinguish these features. 
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In some cases, astrophyslcally Interesting species have transition 
wavelengths which are not precisely known. Such species may be diffi- 
cult to produce in the lab, or may have excited levels from which radia- 
tive transitions are preempted by colllsional de-excitation, further- 
more, the level energies may be difficult to precisely calculate quantum 
mechanically, hence, the most accurate energy level diagrams for such a 
species may be drawn from astrophysical measurements of the transition 
frequencies, . 

2 , 2 * . Inspiration.^ 

There are two categories of high resolution incoherent detection 
techniques In current use in the FIR. The first group, referred to as 
monochromators, samples one spectral resolution element at a time and 
Includes such instruments as grating spectrometers and Fabry-Perot <FP ) 
spectrometers. The second group Includes multiplexing spectrometers 
such as Michel son Interferometers or lamellar grating spectrometers. 
These Instruments sample a number of Independent spectral elements 
simultaneously and recover the power spectrum through Fourier transform 
techniques. We refer to these as FTS (Fourier transform spectro- 

m, • 

•4eal spectrometer would have variable resolution so that one 

resolve lines when one wishes to measure line intensities or 
ovftr-resolve a line when one wishes to obtain Doppler width or Doppler 
shift information. Narrow FIR lines have line widths less than 1 
km/sec, corresponding to a resolution > 3xl0 5 , so one might consider 
construction of an instrument with such resolution. In practice, such 
an instrument is quite difficult to construct. A more reasonable design 
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goal would be R " 10 1 * at which transitions such as the lines due to 
shocked CO In the BN-KL region of the Orion Nebula would be just resol- 
ved and diffuse emission lines from extragalactlc sources would be 
over-resolved# 

The most straightforward technique for achieving high resolution Is 
the construction of a grating spectrometer# The first FIR lines from 
astronomical sources were discovered with such devices# Given back- 
ground limited detection and plenty of work space, such an Instrument Is 
probably more useful than a FP or FTS# A well-designed grating spectro- 
meter may have a throughput approaching 4056, while our own FTS is only 
* 17% efficient and good FP systems pass only 1045% of the radiation 
Incident upon them {Watson 1982). Furthermore, a grating Instrument is 
often readily tunable over a continuous range of wavelengths. For 
example, our medium resolution spectrometer is continuously tunable 
from * 100 to 200 urn In first order, and may be used from - 55 jun - 
200 um if one of the two detectors operates In the second order of the 
grating. Finally, a grating monochromator permits the use of many 
detectors in the focal plane along the direction of the dispersion of 
the grating. It is not unusual for a grating spectrometer to have 10 
detectors in the focal plane leading to a net sensitivity gain of /To 
over a FP Instrument of the same resolution which typically may illumin- 
ate just one detector. However, a simple grating instrument with R - 
lO 4 would require a grating - 75 cm long. Such an instrument would be 
difficult to install on the KAO much less be kept helium cooled. (On 
the Lear Jet such a grating would protrude through the fuselage.) 
Reasonable alternatives for high resolution work in airborne astronomy 


are folded light path Instruments such as Fabry-Perots or interfero- 
metric Instruments whose light path may also be folded and whose reso* 
lutlon achieving drive stage may be kept at ambient temperatures* 

Much of the data presented In this thesis were taken with a Fourier 
transform spectrometer. The actual description of the apparatus awaits 
Section 2.3, but suffice It to say that the FTS system Is a beant*spl1t« 
ting device which acts like a lamellar grating spectrometer with just 
two lamellae* 

This system has several advantages over a comparable FP system plus 
Some disadvantages* The most telling point for such a comparison is the 
source of system noise. Under electronic, or detector noise limited 
conditions, a FTS has a clear advantage over a FP system of the same 
resolving power. This is because the FTS samples many resolution ele- 
ments at once, resulting in the “Fellgett", or multiplexing advantage, 
first pointed out by Fellgett (1951). Individual resolution elements 
are encoded in a Fourier series which is recovered using FT techniques. 
The result (see Section 2.3.2, below) is a gain of /m/8 in the final 
signal /noise ratio over a similar monochromator making the same measure- 
ments, where m * number of measurements (resolution elements). Typical- 
ly we make 32 or 64 measurements, resulting in a multiplexing advantage 
of 2 or 2/2 if the system is electronic-detector noise limited. 

We must stress, however, that the Fellgett advantage does not hold 
under background limited conditions. In fact, at the inception of our 
FTS in present use, our grating spectrometer was clearly pot background 
limited but rather electronic noise limited. Consequently, in the early 
spectra obtained with our interferometer the multiplexing advantage 
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resulted In FT spectra whose sensitivity was two times that which the 
zero path (no interference) signal-to-noise ratios would indicate. 

Since then, with Improvements In detector design and processing, filter- 
ing and electronics we have seen our sensitivity improve by over a 
factor of six. The Fellgett advantage Is now not so evident and it is 
likely that we are approaching background limited performance. 

It should be noted that if the background limit is due to changes 
in emissivity, transmission or turbulence along the light path, it 
should be possible to regain much of the Fellgett advantage through use 
of a reference detector which monitors these changes. If, however, 
photon noise dominates, no such deconvolution is possible and the 
Fellgett advantage vanishes, 

It is also true that the Fellgett advantage is only useful if one 
wishes to sample all of the m points in the FT spectrum. In the pre- 
sence of a continuum one many typically wish to sample eight continuum 
points in addition to the spectral line of interest. Thus, our sampling 
of 32 or 34 points corresponding to 16 or 32 resolution elements re- 
quires that there be two or four spectral lines of interest in the band- 
pass of our instrument. In the FIR this is not an unusual situation, 
the case where there are not two or more interesting lines in our gra- 
ting bandpass of 2.2 pm is the exception rather than the rule (see the 
Appendix). 

In the absence of the Fellgett advantage, the FTS system has 
several other features which make it a desirable system. 

1) The interferometer is a ready extension of our grating spectro- 
meter which acts as an order sorter for the interferometer. We presently 


use our order sorter with 4x4 mm entrance slits which results in a 
resolution R “ 200, 

2) The FTS has advantages of a grating Instrument that a FP system 
doesn’t have. Our FTS system Is fully as tunable as our grating spec- 
trometer } 100-200 pm continuous coverage in first order of our grating. 
FP systems often do not supply full tunabillty, requiring preflight 
selection of which order sorting fixed FP plates will be used. This may 
limit the instrument to about four spectral regions of Interest during a 
flight. A third characteristic our FTS shares with a grating instrument 
is the ability to illuminate multiple detectors. Our grating spectro- 
meter at present contains two detectors working in adjacent wavelength 
regimes. This is a great boon for "search" capabilities, i.e.» If the 
rest frequency of the species in question is unknown, or if we wish to 
sample lines which are closely spaced in the FIR spectrum. 

3) An advantage unique to the FTS is absolute wavelength calibra- 
tion. The wavelength scale of our FT spectra are determined solely by 
the step and burst sizes with Which a mirror is advanced (see Section 
2.3 below). Thus we require no "in-flight" wavelength calibration* If 
the step size is known to relative accuracy A, then the absolute wave- 
length scale of the FT spectrum must also be precise to A. 

4) The FTS shares an advantage with the FP instrument over a gra- 
ting spectrometer in that the spectral resolution may be continuously 
varied without losses in system throughput. Furthermore, only minor 
modifications are required to increase the travel of the drive stage. 

Our initial drive stage had a 5 cm travel corresponding to a maximum 
spectral resolution of .1 cm- 1 . Since then, we have increased the drive 
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stage travel to 16 cm corresponding to a factor of 3.2 Increase in 
spectral resolution. The design modifications were minor, and the new 
drive stage performed quite well on its first flight — confirming the 
existence of the 2 n j/2 (J*3/2 •>• J«l/2) OH emission features in the 
BN«KL region of the Orion Nebula. (Viscuso et al. 1985). 

5) A final calling card for the FTS Instrument Is cost and ease of 
construction. The interferometer was conceived, designed and flown in a 
period of' just seven weeks in March-April 1981. The total cost was 
under $2000. 

2.5. Interferometer., 

2.3,1* Theory* 

Before giving a detailed description of the design and operation of 
the interferometer used for data in this thesis, let us first digress to 
give a simple description of how a Fourier encoding spectrometer works. 
We follow the notation of Harwit and Sloane (1979). 

In the simplest of FTS systems, a collimated coherent beam is split 
evenly across two reflecting mirrors, one of which is movable, call it 
mirror M; the other, being fixed, call it mirror F. Imagine a monochro- 
matic wavefront of frequency, v, and electric field amplitude A incident 
on this system. As mirror M moves a physical displacement, e, the 
optical path length of the beam from M is varied by a distance 2e. 

Thus, the return beam from M lags the return beam from F by a phase 4 - 
4irvs where v now is the wave number. At zero path (e =0), the returned 
amplitude from either mirror is A 0 - A/2. At path difference e, the 
return amplitude from M+F is given by 
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A(e) « A/2 £1 + exp(4ir1ve)] 

Thus, the Intensity at the detector (after the beams have been recom- 
bined) is 

n(e) » A(e)*A(e) ** A 2 /2 Cl + cos(4irve)3 C2.1] 

If the Incident wavefront is the superposition of many wavenum- 
bers, v, each with its own intensity, then the total intensity in the 
frequency Interval v + v+dv may be written i|>(v)dv so that 

n(e) a X/2 / tp (v ) CX * co$(4itve)] dv 
o 

a 1/2 + 1/2 / ip(v ) cos(4wvs) dv [2.2] 

o 

where s / tp(v ) dv is the intensity at zero path. The 
o 

interferogram thus consists of a modulated part (the integral) super- 
posed on a fixed O.C. offset. The physically interesting part is the 
Integral: 

n'(e) » 1/2 J Ip(v ) cos(4rrve) dv 
0 

Since i(i (v) and cosine are symmetric functions, we may write: 
n'(c) .» 1/4 / i(/(v) e 47rl ’ VE dv 

•-ctf 

where we have replaced cps(4irvs) by e^ ive as the anti -symmetric 
sine term will vanish. Thus defining x s 4 tiv we have: 

ti'(e) - I/16 tt /* ip (x/4ir ) e^ Xc dx 


We recognize the Integral as the Fourier transform of *Kv) s ¥(47re), 
thus: 

n ' (e) « 4tr/167r ¥(4ire) or, 

V (47re) *» 4n*(e) 

Thus the signal at the detector, n'(e), Is proportional to the 
Fourier transform of the Incident waveform tf(v) 

X/4ir / W T(4ire) e" 4irjsv % de 
* 4/ n'(e) e" 4 ” 1ev de 
* 2ir<Kv) a Iffir / n * (e ) s m4it ^ sv de 

*40 

or again appealing to symmetry: 

ip(v) a 16 | n'(e) cos(4irev) de 
0 

To precisely recover tfi(v), we must sample n(e) continuously from 
zero-path to Infinity and perform the cosine transform of n^e). In 
practice, of course, it Is Impossible to sample the entire Interfero- 
gram. Instead, we sample n(e) at n discrete, evenly spaced points a 
distance A apart along our total travel of (n-l)A, The step size. A, is 
chosen such that A a « v m j n )], where v^x is the maximum 

frequency to be sampled and v ml * n is the minimum sampled frequency (see 
Section 2.3,3 for comments on choosing step size and its relation to the 
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Nyqulst sampling rate). From this point on, we drop the superscript 
from n'(e). 

The Interferogram thus sampled may be discretely transformed to 
form an estimate $j(v) of the true spectrum 

a n m 1 

$i(v) ** 16 I n(rA)Acos(4tfrAv) 
r«{J 

We wish to evaluate $(v) at the n discrete frequencies v s given by 
v s K v min * for' s ~ 0,l»».. t n M l and 6 s ( v max** v m1n^^ nM ^* 

The true estimator of if»(v) in the frequency range v m j n ^ v s s v max 
is given by; 

$(v) ** $i(v)6 « $i(v)/4(n~l)A or, 

$(v) s 4/(n*l) l n(rA) cos(4irrAv) [2,3] 

r=0 

We formally rewrite this as the matrix equation; 


¥ « A n 


where £ and n are the column vectors defined by? 



*( v min) 


n(0) 


$( v min +$ ) 


n(A) 

Y a 

f 

« 

n * 

. P . 

• 


P 

$( v max) 

»+ •*! 


nC(n4)a] 

ft 


and A Is an nxn matrix whose elements a sr are defined by 
a$r « 4/(n«l) cosC%rA(v m ( n +s6}] 

for s,r » 

2*3,2. The Fellgett Advantage.. 

The factor of 4/{n*4) in equation 2.3 above is very Important 
This is the term which leads to the Fellgett advantage. 

The Fellgett advantage is realized if the following criteria 
apply: 

1) The dominant source of noise is detector-electronic 
noise, 1.e M the noise is Independent of photon flux. 
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2) The noise of each measurement:, ej, has zero mean? 

<e-j> * 0 (where < > indicates a time average 
expectation value)* 

3) The noise between various measurements 
n(e|), n‘(ej) is uncorrelated, i.e., <ejej> * 0. 

4) The noise has a variance defined by <(ej} 2 > » o 2 . 

Under these circumstances, the expectation value of the estimator f 
is given by: 


<£> « t(v) 


and we may write: 

* - 1 ) 1 « A e + * a s>0 e 0 +...+ a s>n 4 e n „i 

where e is the noise column vector. 

We define by s Vg the variance of the spectral estimator for 
frequency v$» T(v s ) 

e Vs s <E$(V$) - f(v s )3 2 > « o 2 (a| j0 + ... a| >n «l) 


the average mean square error, e, is given by: 


e 5 (1/n) l e Vs « cr 2 /n (aJ j0 + af |X + a^l,n-l) 

« u 2 /n Trace (A* A* 1 * ) 

B 16a 2 /[n(n~l) 2 ] l £ cos 2 C4irrA(v m i n *sS)3 

r=0 s«0 

~ 8na 2 /(n**l) 2 + 8a 2 /Cn(n-l) 2 3 l X cos[8irr4(v m j n +sS)] 


" So^/n for large n. 

Thus, If we measure a variance o 2 for each n(e), then we will measure a 
variance In the transform spectrum which Is smaller by a factor of 
(8/n) 1 ^ 2 . This Fellgett advantage Is quite powerful when a detec- 
tion system Is detector-electronic noise limited* For a discussion of 
the relative merits of multiplexing Instruments, see Harwit and Sloans 
(1979). 

Prior to the Introduction of J-FETs and better detector processing 
techniques, our FTS system was clearly detector-electronic noise limi- 
ted. Figure 2.1 shows the Fourier transform spectrum representing the 
first astronomical detection of the d » 16*15 rotational transition of 
highly excited CO (Stacey et al. 1983). This spectrum Is the average of 
three spectra taken in January 1982, while observing the BN-KL region of 
the Orion nebula. For each spectrum, we sampled the interferogram 33 
times and Integrated 10 seconds at each point. Thus, we expect to see a 
factor of two (« (n/8) L/f2 ) Increase In sensitivity In the Fourier 
transform domain relative to raw data* At zero path we measure a system 
NEP, 1 Including all losses, of ~ 3.6*10" 13 w The grating 

bandpass Is roughly Gaussian In shape for an extended source such as 
BN-KL. The dashed line represents a least squares Gaussian fit to the 
data, Ignoring the telluric water vapor absorption features at 161,90 pm 


1 Noise Equivalent Power (NEP) is the power incident on the telescope 
which produces a signaT/noise ratio of one, in one second of integra- 
tion. 
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Figure 2,1 Fourier transform spectrum of the BN«KL region of the I 

I 

• • . • • . . • - I 

Orton Nebula showing the first astronomical measurement of the j 

J « 16*15 (162.81pm) rotational transition of highly excited CO j 

(Stacey et al. 1983), The spectrum has been fit with a least squares 

Gaussian centered at 162,55pm which accurately mimics our grating 

profile. The CO feature and nearby H 2 0 features at 161.90pm and 

163.53pm have been ignored in the fitting process. 
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and 163.53 pm, and the 162.81 urn CO feature. If one measures the 
deviations of the data from this best fit, one calculates a system NEP 
* 1.9»10" 13 W Hz" 1 ^ 2 "•* a factor of two better than the NEP of the 
Interferogram — as predicted. 

Recent advances In electronics and detector processing have 
Increased our sensitivity so that our Interferometer - spectrometer 
system is now typically characterized by system NEP l s * 5 h 10“ 14 W 
Hz'* 1 ^ at zero path. Recent spectra show that the Fellgett advan- 
tags is much smaller than predicted -« resulting In less than 20S» 
Improvement In the Fourier transform domain. We take this as good 
indication that our system is no longer detector-electronic noise 
dominated and approaches the sensitivity limit imposed by random 
background fluctuations. 

2,3,3* Choosing Burst Size: The Nyquist Rate, 

If we wish to sample all the frequencies, v, of <Kv) from. the v « 

0, D.C. verm to some maximum frequency, v mx , we must sample rj(e) at 
the so-called Nyqui st rate to avoid aliasing. 2 The Ityquist, or samp- 
ling theorem (see, for example, Papoulis 1977), simply states that to 
obtain an unaliased spectrum in the frequency domain, we must sample the 
Interferogram at a rate equal to twice the highest unaliased frequency. 


2 Aliasing refers to uncertainty concerning the spectral order of the 
interferogram, or transfer from one order to another. 
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For the Interferogram, this means we must sample n(e) every l/4v max 
cm. (Remember, every physical displacement e corresponds to an optical 
path length Increment 2e.) For a drive stage 5 cm long at a frequency 
v max ’ 67 cm** 1 ( 15 ° Mm) this means we must sample n * (e ) over 1000 
times for just one spectrum. Clearly this would be prohibitive in terms 
of both integration time and the computer time needed to perform a 1000 
point Fourier transform. To both circumvent this problem, and especial- 
ly to reduce the background flux incident on the detectors, we employ 
our high resolution spectrometer (see Section 2,4, below) to act as an 
order sorter and greatly reduce the number of frequencies to be sampled. 
At 150 pm, our two-detector spectrometer samples a total bandpass of 
roughly 2 pm corresponding to Av spec t - .9 cm* 4 . Thus we need only 
sample fast enough such that we sample all frequencies between v m ^ n 
and v max where v min is determined by v min < Av spect $ v max . 

This may be ensured by choosing the burst size, a such that a = 
l/[4(v max «v m *i n)3* 

In practice, we choose a by the following convention. First, 

choose the central wavelength, v P » of the Fourier spectrum. Typically 

this wavelength is chosen such that it falls halfway between the peak 

spectral response of the two channels in the spectrometer. The central 

wavelength is determined by v c = (m+l/H)/4a, where m * interferometer 

order, usually - 30-40 for our 5 cm travel stage. Intuitively, m 

denotes the number of interfering wavefronts skipped with every burst. 

The infrared frequency v c is aliased through this technique to a much 

lower frequency \)' c = v c /m. The factor of 1/2 above, then, ensures 
* 

we sample each aliased frequency twice per wavefront. Now A Is 
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physically constrained by the drive stage travel such that if we 
choose to use n bursts to sample our spectrum, we must have nA £ 5 cm. 
Typically n « 32, so that a - 1500 ym. The free spectral range of the 
resultant spectrum is given by the difference between adjacent Inter* 
ferometric orders or, Av^ e v c z [4A(l/m * l/(m*l)V Having 
chosen the step size, 4, one may assign frequencies to the Fourier 
transform spectrum using the reverse of the above transformations. 


2»3.4. Functional Description* 


Our interferometer was first described in a paper by Harwit et al. 
(1981). The interferometer acts like an all-reflecting lamellar grating 
spectrometer with just two lamellae (see Figure 2*2). The infrared 
radiation from the KAO's f/17 beam is reflected off a dichroic and 
brought to a focus at plane mirror A. The visible light passes through 
the dichroic and is imaged within the offset guidar reticles. (Designed 
by Prof. James R. Houck.) At mirror A is a 4 mm wide entrance slit 
which ensures our sampling within the primary diffraction lobe of the 
telescope (~ 6 mm from 1 st null to 1 st null at 150 um) and thus 
sampling a coherent beam. Mirror A diverts the diverging beam to plane 
mirror B which reflects the light into C. C is a 45 cm focal length 
off-axis paraboloid which collimates the f/17 beam into a beam - 2.5 
cm in diameter. The collimated beam then strikes plane mirror D which 
is aligned so that the circular beam is split evenly across mirrors E 
and F» Mirror E is the plane mirror whose movement relative to plane 
mirror F causes the interference. These mirrors have been aligned so 



Figure 2.2. Our lamellar grating interferometer top view (top) and side 
view (bottom). The dashed line shows the light path through the system 
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that at zero path there Is just a fraction of a millimeter between the 
two, thus minimizing the effects due to diffraction. As E moves, the 
path length, 2s, of the return beam Is varied. Figure 2.2 shows that 
the path of the collimated beam is parallel to the direction of movement 
of the drive stage and above mirrors C and 6. Thus, there is no 
blockage of the light from F as E moves to full path difference. The 
beams from mirrors E and F are recombined at P before striking G, n 
second 45 cm focal length paraboloid which refocuses the beam to f/17 
and images it on the entrance slit of the spectrometer. 

Note that when mirrors E and F are coplanar, the spectrometer 
receives exactly (minus losses due to the seven reflections - 30ft) the 
same radiation it would receive from the telescope in the absence of the 
interferometer. 

A 5 cm travel drive stage was used for all of the interferometer 
data presented in this thesis. The drive stage is commercially 
available (Aerotech, Inc., Allison Park, PA 15101) with an accuracy 
better than 1 part in 10 3 and with somewhat better reproducibility. The 
stage is driven by a stepper motor and geared in such a way that one 
step of the stepper motor corresponds to a 2 um translation of the drive 
stage. We were able to accurately calibrate the step size through our 
detections of the 0 « 17+16 and 0 * 16+15 rotational transitions of 
highly excited CO found in the BN-KL region of the Orion nebulae. The 
wavelengths of these transitions are known to very high accuracy (Todd 
et al, 1976). These detections determined that our step size was actu- 
ally 1.998 ym — within the advertised accuracy of the drive stage. 


The spectral resolution of our drive stage is determined by the 
longest optical path length difference introduced into our interfero- 
grant* For a 5 cm travel* the optical path difference is 10 cm corre- 
sponding to a spectral resolution of .1 cm" 1 * At 157 urn, this corre- 
sponds to a resolution R s X/AX - 650. 

This resolution estimate is based on the standard Rayleigh criter- 
ion. That is, we will resolve two lines of equal intensity as indepen- 
dent spectral elements if they are separated by av - .1 cm" 1 . The 
full width half maximum (FWHM) of any unresolved spectral line may be 
made a factor of two narrower than this through proper choice of burst 
size. We define our spectral definition, D, as 0 s x/axfmhm * 2R. 

Thus, we may measure line widths and wavelengths a factor of two better 
than the resolution* The definition of our 5 cm travel stage is D - 
1300, or Av * 230 km/sec. 

To fully achieve this resolution, one must sample the interferogram 
at both positive and negative path lengths. However, if one a priori 
know: that the interferogram has begun with a sample at zero path, then 
one may mathematically reflect the data about zero path, as we are 
assured by symmetry that signals- measured at a positive path difference, 
e, must be the Same (within the noise) as those measured at negative 
path difference, ~e. In practice, determination of zero path has been a 
straightforward process. One first eyeballs the movable mirror, E, to 
within a fraction of a millimeter of zero path on the negative side. 
Then, with the .grating spectrometer set at a suitable wavelength, 
typically 118 um, one single steps the movable mirror through zero path 
producing a sinusoidal fringe pattern. If one then repeats the process 
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at two additional widely separated wavelengths, say 140 ym and 163 ym, 
avoiding rational relations between wavelengths, one may then unambig- 
uously determine which fringe peak Is common to all wavelengths* This 
fringe peak must then be zero path and zero mirror-position Is then 
redefined as a convenient number, usually as step 1000, and used as the 
reference point throughout the flight* 

0 

2,3*5. Design limitations.. 

The efficiency of our interferometer is strictly governed by the 

modulation depth, i.e», the fringe visibility near zero path. As can be 

seen from equation 2.2* for a sufficiently narrow bandpass, as one steps 

the movable mirror about zero path, the two beams will interfere to form 

a sinusoidally varying interferogram of maximum power P « and minimum 
■* 

power 0. This is the ideal system,' let us denote its efficiency by e » 
1. In reality, various diffraction effects lower the coherence of the 
recombining beams. This results in both a lowering of the peak inter- 
ference signal, P, and an increase in the valley signal, V, The effi- 
ciency, e, may then be defined as e s (P-VJ/P, We can see that if there 

is no modulation, i*e,, no interference of beams, P = V and e = 0. For 

* 

coherent lab sources, we have obtained fringe visibility ratios P: V as 
high as 6:1 resulting in an 83% efficient interferogram. Typical In- 
flight fringe visibility is - 4:1 resulting in a 76% efficient scan. 

It has been found that the fringe visibility is extremely sensitive to 
the gap between mirrors E and F at zero path. Fringe visibilities as 
high as 6:1 with a gap Of .1 mm are degraded to - 2:1 if the gap is 
allowed to grow as large as ,5 mm. 
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An unfortunate consequence of the original design for the Inter- 
ferometer Is that mirror E had to be aligned such that Its reflecting 
surface was not perpendicular to Its direction of travel. Thus, the 
collimated beam walks across mirror G as E travels and consequently also 
moves across the secondary of the spectrometer. Portions of the beam 
eventually walk off the secondary resulting In a loss of light from E at 
long path lengths, hence a loss of both modulation depth and absolute 
signal level* An interferogram of a strong continuum source (BN-KL in 
Orion) then has the characteristic “droopy shoulder" appearance of 
Figure 2.3a. The Fourier transform of such an Interferogram Invariably 
shows a huge first cosine term. The data may be mathematically correc- 
ted by fitting a least squares line to the interferogram and removing 
this baseline from the data (Figure 2,3b). Recently, this problem has 
been alleviated entirely through introduction of corner mirrors. The 
return beams from corner mirrors E and F are then both parallel to the 
direction of travel of E resulting in the elimination of walk-off 
problems. 

2 Am Spect rometec*- 

2.4. 1» Layout and Functional Description. 

The spectrometer used for the observations presented in this thesis 
is a fully liquid helium cooled grating spectrometer operated in a 
Littrow mode. Essentially, it is an improved version of the spectro- 
meter originally described by Houck and Ward (1979) and later improved 
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Figure 2.3 a&b. Apodization of an interferogram. Figure 2.3a shows 
the average of four i nterferograms taken of the BN- KL region of the 
Orion Nebula of September 23, 1981. Zero path is the high point on 
the left and the point 33 on the far right is data taken at largest 
travel (5cm) of our drive stage. Error bars are one standard devi- 
ation from the mean. Figure 2.3b shows the same interferogram 
following droopy shoulder removal and apodization with a triangle 
(Bartlett) window. The D.C. offset has been subtracted out. 


Wavelength (microns) 


Figure 2.4. Fourier transform spectrum of the BN-KL region of 
Orion. This spectrum is the Fourier transform of Figure 2.3b } 
follov/ing reflection about zero path. The fitted Gaussian 
represents the grating instrumental profile. The feature at 
153.3pm is the first astronomical detection of the J « 17+16 
rotational transition of shocked 12 C0 (Stacey et al. 1982). 
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and reconfigured for the f/17 beam of the KAO by Melnlck (1981). As can 
be seen In Figure 2.5 the collimator Is a Cassegrain system. Light 
entering the spectrometer at A passes through the far Infrared filter at 
B before being deflected to the secondary by a 45° plane mirror at C. 
Just In front of C Is the entrance slit. C Is set slightly off-axis so 
that the return diffracted beam misses the mirror. The secondary Illu- 
minates the 3" primary, 0, which then sends the collimated beam to the 
diffraction grating at E. The diffracted beam then is de-col 11 mated and 
exits the collimator at F. The collimator bolts directly onto the 
aHe-cooled wort surface, which Is the bottom of the AHe can. The beam 
at F Is deflected by two 45° mirrors displaced - 4 mm along the 
direction of dispersion on the work surface. These mirrors send the 
beam onto the exit slits located at the entrance to the light cones of 
the two detectors. All of the structural parts of the collimator are 
constructed of stress-relieved aluminum. 


2.4.2. Optics.. 

The effective focal ratio of the collimator may be adjusted from 

f/6.5 (to match the Lear telescope) to f/17 (to match the KAO telescope) 

by changing optics. The f/6.5 optics were made of low thermal expansion 

Pyrex glass by 3-B Optical Co. (Gibsonia, PA). Both the concave primary 

(radius of curvature, R c - 38,9 cm) and the convex secondary (R c “ 

24,15 cm) have spherical figures as the diffraction lobe of both the 

Lear and KAO telescopes is large enough at 150 pm to make spherical 
* 

aberration negligible. Glass mirrors have two inherent problems. They 
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have a large thermal mass and time constant and shrink differentially 
relative to the aluminum collimator walls. Thus* the focal length of 
the Instrument must be adjusted at room temperature so that the focus of 
the Instrument shifts the correct amount as the collimator cools to 
liquid helium temperatures. This differential contraction results in a 
shift in the focus - 8 mm. The advent of high surface quality pure 
aluminum mirrors has eliminated both of these problems. Our f/17 optics 
are all aluminum, primary R c * 27.06 cm, secondary R c « 5.76 cmj 
machined by a diamond tool at Aeroresearch Corporation, Great Neck, NY. 
Essentially all of the f/17 collimator is then aluminum, which means the 
unit shrinks uniformly and therefore the focus does not shift as the 
collimator cools. The aluminum mirrors have the further advantage that 
they may be directly heat-strapped to the collimator walls through 
copper braiding. These straps are attached with epoxy or screws. 

The all -aluminum grating (Bausch and Lomb, Rochester, NY) has been 
blazed at 26.75° and ruled at 6 11 nes/ntn for first order work at 150 
pm. The grating also has been heat-strapped to the collimator walls. A 
stepper motor, external to the helium dewar positions the grating. The 
stepper motor's motion is geared to a fiberglass rod which feeds through 
an Q-ring in the dewar bottom. The fiberglass rod minimizes the heat 
load into the dewar. The grating is coupled to the motor through a set 
of bevel gears and a worm-wheel assembly. The stepper motor steps 
through 7000 steps to provide a 30° rotation of the grating. This 
corresponds to grating tunability from ** 100 pm < X < 200 pm, each 
step corresponding to - .019 pm at 150 pm. The two channels are 
separated in wavelength space by one resolution element or - .7 pm in 
the f/17 configuration. 


Wavelength calibration of the grating spectrometer Is achieved 
through matching observed water Yapor features In the spectra of 
laboratory continuum sources with the known wavelengths of these 
features {Traub and Stler 1976). Our stressed OesGa detectors are 
sensitive from 100 urn $ \ $ 195 pm and we generally Identify a dozen 
dozen or so water. vapor features within this regime with which we call" 
brate wavelength as a function of grating step. Theoretically one 
should use the grating equation to make the calibrations which, for 
Ltttrow systems reads: 

mx/d * 2 sin(3) 

where m « order, d » groove spacing, and $ « tilting angle of the 
grating with respect to the normal* Converting tilt angle to grating 
step Is somewhat tedious, and It Is found that In practice, a simple 
linear interpolation procedure between water vapor features Is suffl « 
dent. The accuracy of our grating calibrations, as verified with the 
Interferometer, is typically better than #2 pnt or 1/4 resolution 
element. 

For Lear data, these grating calibrations, as verified by In-flight 
spectra of calibrators, are used for the actual wavelength calibrations 
of the final spectrum (see section 2,6,1 below). As the interferometer 
is self ^calibrating (see section 2.3.3), we generally take interfero- 
grams of laboratory continuum sources before take off on KAO flights to 
verify our grating calibrations and make any necessary adjustments. We 
thus ensure that the grating bandpass Is centered on the expected line 
position and well within the free spectral range of the interferometer. 
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As long as the stepper motor colls are energized, the spectrometer 
Is reproducible from scan to scan, to within a step. Between flights, 
the grating Is re«zeroed (set to the short wavelength mechanical stop) 
by hand, a process which introduces a flight to flight reproducibility 
error of a few steps, or • .04 pm. As the spectrometer resolution Is 
- .7 pm, this error is negligible* 

Diffraction determines the smallest logical slit size to use for 
astronomical observations, hence maximum spectrometer resolution. The 
30 cm f/6.5 telescope aboard the Lear Jet facility is. diffraction 
limited at 150 pm to an Airy disk - 2.2 mm in diameter (null to 
null). This then Is the optimal slit width along the direction of dis* 
perslon. For a 2.2 mm slit, we have measured system resolutions for 
point sources of R s x/tx - 200. The 91 cm f/17 telescope aboard the 
KAO has a diffraction lobe * 6.2 mm in diameter. In flight, we 
generally have used a 4 mm wide slit, which results in a small ( * 1%) 
loss in luminosity. The spectrometer in this configuration also shows a 
resolution of about 200, 

2.4.3. Dewar, 

The spectrometer and helium can are contained within a 5" diameter 
cylindrical dewar. The helium cooled spectrometer and helium can are 
shielded from external radiation with an internal heat shield which fits 
neatly over the spectrometer helium can assembly and is attached to the 
top of the helium can* The heat shield is wrapped with ~ 20 layers 
of superinsulation (aluminized mylar) and maintains a temperature 


somewhat lower than 5Q°K, The hold time of the 3.5 liter helium can Is 
* 28 hours unpumped and * 20 hours when the helium coolant Is pumped 
down to 2°K. The dewar Is sidelooking and the vacuum seal Is preserved 
through the light port by a *010" thick piece of clear polyethylene. 

The polyethylene Is rugged and transmits over 90% of the FIR radiation. 

2*4*4. Filters* 


As our Ge:Ga detectors are sensitive to near infrared radiation, 
especially - 2 pm where the germanium is an intrinsic photoconductor, 
the collimator heat shield assembly contains a series of filters to 
block unwanted radiation. The visible and near infrared are quits 
effectively blocked by a 5 mil thick piece of black polyethylene; i.e., 
polyethylene which has been Impregnated with carbon. Since the black 
polyethylene scatters as well as absorbs the short wavelength radiation 
(McCarthy 1967), the filter is placed on the inside of the heat shield, 
well away from the detectors. Virtually no light Is transmitted from 
2,5 - *7 pm, where the black polyethylene slowly begins to pass, 
reaching a transmission of ~ 8% at ~ 50 pm. At 150 pm, the poly- 
ethylene transmits - 70%. 

To block unwanted radiation longward of - 32 pm, we employ two 
helium cooled reststrahlen filters, The term reststrahlen (residual 
rays) is applied to certain salt crystals which display a high degree of 
reflection (£ 90%) at wavelengths corresponding to the ionic vibrational 
energy bands in the crystal lattice. In the entrance cone of the spec- 
trometer we have a 1 mm thick piece of BaF 2 , which effectively blocks 
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from 18 - 58 pm (Armstrong and Low 1973). The BaF 2 has been anti reflec- 
tion overcoated (IR Labs, Tucson, AZ) for work at 157 ym at which 
wavelength It transmits - 90%. At the exit slits, directly in front 
of the detectors we have Installed a second reststrahlen filter, a .5 mm 
thick piece of AgCl. AgCl passes - 80% at 157 ym and effectively 
blocks 42 ym $ X < 100 ym (Armstrong and Low 1973). An .alternative 
far-infrared reststrahlen filter would be a combination of CsBr and KBr 
filters In place of the AgCl (Watson 1982), These filters block 52 $ X 
$ 130 ym and also transmit - 80% at 157 ym, However, these filters 
are much more delicate than the AgCl. They both are highly hygroscopic 
and KBr 1$ quite susceptible to the thermal shock of cooling down the 
spectrometer* AgCl has the distinct disadvantage of possessing a deep 
absorption feature near 152 ym which the CsBr an<2 KBr do not possess 
(Figure 2,6). Furthermore, AgCl filters are sensitive to UV radiation 
and they slowly oxidize on extended contact with the brass light cones. 
The former problem is alleviated through care to minimize UV exposure, 
while the latter is eliminated by isolating the filters from any metal- 
lic contact, with Teflon shim stock. For a description of the manufac- 
ture and implementation of reststrahlen filters, see Watson (1982). 

The total transmission of our spectrometer, including losses due to 
imperfect reflection off the collimator mirrors (10%), grating efficien- 
cy, and filter transmission efficiency is roughly 20%. 

Stray light can be a major problem in a Littrow mode spectrometer, 
as baffling is difficult when the entrance beam is so close to the exit 
beam* The "channel 1" detector is particularly susceptible to stray 
light reflected directly back from the secondary, as this detector Is on 



Figure 2.6. Laboratory grating calibration spectrum of liquid N 2 
cooled glass beads through two different sets of reststrahlen filters. 
Most of the absorption features in these continuum spectra are due to 
atmospheric water vapor absorption in the 50cm light path from the 
spectrometer to the glass beads. The first panel shows the spectrum 
obtained when CsBr and KBr are used as cut on filters to pass beyond 
-135pm.- The second spectrum is taken with AgCl which passes longward 
of -100pm. Note the-AgCl absorption feature near 152pm. No effort 
has been made to deconvolve relative detector or filter response, The 
varying depths of H 2 0 features from the first to the second spectrum 
are due to changes in laboratory relative humidity from day to day. 
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the optical axis. This problem has been minimized through placing a 
cone of black carborundum sandpaper in the center spot of the secon- 
dary. Baffles have been placed wherever possible in the spectrometer, 
Including an entrance baffle, a baffle at the primary and one about the 
secondary. The entire Inside of the spectrometer has been painted black 
with Nextel brand black velvet spray paint mixed with carborundum. This 
combination reflects roughly 10 % of the FIR. Residual stray light in 
our system is of the order 5-40%, 


2.4.5. Electronics- 


Changes in resistance of our Ge:Ga photoconductors are converted to 
changes in current through a standard transimpedance amplifier (TIA) as 
shown in Figure 2.7. The primary advantage of a TIA is that it main- 
tains a constant voltage drop, across the detector which ensures 
linearity of system response. The electronics may be broken into three 
temperature stages: the detector resistance,- Rp, and load resistance, 
at 2°Kj the preamplifying J-FETs at 70°Kj and the room temperature 
Operational Amplifier (OpAmp). The matched J-FETs (gain - 1) first 
stage of the amplifier is kept as close as possible to the detector and 
load resistor. This minimizes the length of the high impedance wire to 
the input stage, gate G 1} and hence, minimizes microphonics. 

The gate G 2 is kept at ground, so the action of the OpAmp is to 
keep the gate Gj, at ground. In operation, bias applied across the 
detector, Vp, acts to move free carriers through the detector resul- 
ting in a current flow, Ip. To keep G ls at virtual ground then, the 
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OpAmp supplies a yoltage V oll t across the load resistor such that 
v out s m ^d^L * mV B r I./ r D* The output voltage from the OpAmp 
Is the measured signal* Changes In detector resistance thus result in 
changes in V ou t* One may define a gain from the TIA as G - 
"■AVout/AVg * -Ri/ARp. 

Sources of noise In the TIA circuit are the Johnson noise due to 
the load resistor) detector-electronic noise, and microphonic pickup. 
Johnson noise* due to random thermal motions of electrons is given by 

AVj a (4kTR^Af ) l/2 volts 

where k “ Boltzmann's constant, T * temperature (°K) and Af is the 
electrical bandwidth. Taking T s 2°K, and Af ~ 1 Hz, we have aVj « 
6x10~ 7 (Rl/ 10 9 )^ 2 volts. As Johnson noise scales as \ l ^ z 
and the gain, G, {hence the signal) scales as Rl, a suitably large 
gain will eliminate Johnson noise* In practice, a gain of 1 is suffi- 
cient to make Johnson noise an order of magnitude smaller than back- 
ground noise. Our stressed detectors have a resistance * 10 9 ft under 
operating background conditions, hence we use a TO 9 n Victoreen load 
resistor. At zero bias, we measure Johnson noise. 

The switch-over from MOSFET preamplifiers to J-FETs and improve- 
ments in detector manufacturing techniques has effectively eliminated 
detector electronic noise as the primary source of noise for our sys- 
tem. As work space is severely limited in our dewar, and MOSFETs run at 
2°K right on the work surface, the changeover to J-FETs was a slow and 
painful process, Paul Viscuso was instrumental in implementing this 
change-over. Much of the data taken for this thesis were taken with 
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MOSFETs which effectively limited our minimum NEP to - 2*10” 13 W 
Hz" 1 / 2 . J-FETs need to operate at temperatures above 50% and' as our 
dewar has no available &N 2 surface, the d-FET must be thermally isolated 
from the &He surface and heated externally. At George Gull's sugges- 
tlon, we stand the d«FETs off the work surface with a 1" long "L"»shaped 
piece of thin walled 1/4" diameter fiberglass tubing# This required 
some modification of the spectrometer * helium can interface. The dual 
d-FETs are fit Into the end of the "L" section with a brass ring and 
heated by a 1.1 kfl Allen-Bradley resistor. Thermal contact between the 
heater resistor and the d-FET is assured as the wire running between 
the drains of the two d*P£Ts is wrapped several times around the resis- 
tor and insulated in a "cocoon" of rubber cement. Care must be taken at 
all times not to thermally ground the d-FETs# The fiberglass tubes have 
been notched to minimize heat conductivity and only fine (2 mil) Formvar 
wire is used to make electrical contacts# The wires to the d-FETs and 
heaters are tied together to minimize microphonics and wound several 
times about the fiberglass tube to lengthen the thermal path. The wire 
bundle is tacked in several places along its travel with rubber cement 
and pulled tight to minimize microphonics# The high impedance lead 
between D]_ and the detector is the most susceptible to microphonics. 
Consequently, thicker wire (3 mil) is used for this lead and care is 
taken to minimize its length. 

Whenever possible, all wires are isolated from the work surface 
with Kirn-Wipes paper tissue and tacked with rubber cement. The 
Victoreen load resistor is mounted on a Kim-Wipe with rubber cement as 
close as possible to the detector and J-FETs. The entire system, as 
such, is usually insensitive to mi crophonic pickup. 
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2.4«g« Stressed Detectors. 

The best detectors In current use for work In the FIR are extrinsic 
photoconductors made by low-level doping of germanium with various 
Impurities from columns IIIA and VA of the periodic table* Gallium 
doped germanium (Ge:Ga) was the material used for the first astronomical 
detection of a FIR fine structure line, the 88 pm [OIII] line detected 
by Ward et al. in 1975* 

The spectra presented in the present work were obtained with 
stressed GeiGa photoconductors. High purity germanium Is doped with low 
levels {" 10 llf enf 3 ) of gallium during crystal growth. Germanium has 
four valence electrons and thus forms covalent bonds with four 
neighboring germanium atoms by sharing an electron. Gallium has three 
valence electrons, so at an impurity site there is a deficiency of one 
electron per atom in forming the crystalline bonds. The result is a 
hole which behaves like a positive charge loosely bound to the 
negatively charged gallium atom. In effect, the impurity introduces 
vacant discrete levels, often termed impurity acceptor levels, slightly 
above the valence band. Photons of sufficient energy may then excite 
electrons into the impurity levels leaving free holes in the top of the 
valence bands. In the presence of a bias field, the hole may circulate 
through the electronic circuit allowing current flow until the hole 
eventually recombines at an impurity site. The photoconductive gain, 
Gp C , is defined as the number of passes a hole makes through the 
defector before recombination. A second important quantity in our 
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photoconductor is the quantum efficiency, t)» defined as the ratio of the 

number of Incident photons which free a charge carrier, to the total 

number of Incident photons, In good detector material, impurity levels 

are kept quite low to avoid impurity banding effects and maximize free 

carrier lifetimes, i»e,, photoconductive gain. However, lowering the 

dopant concentrations effectively increases the mean free path for an 
* | 
Impurity site ionizing photon, hence lowering the quantum efficiency, A -I 

solution is to increase the detector's physical size. For typical 

dopant levels the mean free path of an Impurity-site ionizing photon is 

- 5,0 mm, hence the detectors would need to be roughly 5 mm cubes. 

Unfortunately, this also Increases the detectors cross section for \ 

j 

ionizing radiation such as cosmic rays. At aircraft altitudes this can , 

be a serious problem (see Houck and Briotta 1982), so detector dimen- 
sions are generally kept small (1x1x3 mm), The detector then must be 1 

placed in an Integrating chamber to keep the absorption probability 
high. 

Germanium doped with “2X10 11 * cm** 3 of gallium produces acceptor 
Impurity levels " 10** 2 eV above the valence band. Hence, these detec- 
tors are sensitive to FIR shortward of \ «* 130 pm. The application of 
a uniaxial stress along the P-00] crystallographic axis splits the four :■ 

fold degeneracy at the valence edge in the germanium into two doubly j 

degenerate levels (Price 1961, Hall 1962), The splitting is propor- 
tional, to .the strain amplitude, so that a stress of ~ 60 kgf/mm 2 will 
bring the upper valence level to within ,006 eV of the acceptor level. 

The result is photodetective response out to - 205 pm (Kazanskii, 

Richards and Haller 1977; Haller, Hueschen and Richards 1979). 


Stressed GeiGa detectors must be kept at very low temperatures to 
prevent phonon excitation of holes Into the acceptor Impurity levels. 
Typical operating temperatures are T - 2.5°K. 

Experience has shown that most of the stress on a stressed detector 
is due to the differential contraction between the detector housing and 
the germanium crystal. Thus, the housing should be constructed of a 
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material with a large coefficient of thermal expansion. Aluminum imme- 
diately comes to mind. However, housings made of A1 have shown a high 
degree of susceptibility to cold flow, and It may take many thermal 
cyclings for such a housing to reach a state of equilibrium {Watson 
1982), Hard brass appears to be a better candidate for housing material 
as It shrinks nearly as much as aluminum with decreasing temperature, 
but Is less susceptible to cold flow. 

Figure 2,8 shows our design for the detector housing. To minimize 
cold flow and relaxation of the housing, care has been taken to maximize 
the cross section of the brass housing through its weakest point. The 
detector is compressed between two polished stainless steel plugs by a 
1/4-28 set screw « ball bearing mechanism. The ball bearing serves to 
decouple the torque due to the set screw. The Integrating cavity Is 
highly polished and is a snug fit for the stainless plugs Insuring that 
the faces of the plugs are perpendicular to the cavity walls and, hence 
the faces of the GeiGa crystal. Stainless steel was chosen as the plug 
material since it is both very hard and has a small thermal expansion 
coefficient. This then maximizes the differential contraction between 
the housing walls and the steel plugs, hence, maximizing the mechanical 
stress on the cooled detector. 
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Figure 2.8. Stressed detector housing and cutav/ay of finished 
detector, With the light cone removed, detectors are inserted through 
the cavity behind the light cone. All parts are made of brass save the 
stainless steel plugs at either end of the detector and the steel ball 
bearing and set screw. 
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Detector material was purchased from Eagle-Picher in the form of a 
2 3/4 inch diameter wafer cut 3 mm thick along the [100] axis. The 
germanium had been doped with 2»lxl0 14 cm" 3 of gallium with less than 
IQ 13 cm" 3 of residual impurities. Ohmic contacts were created as 
described by Haller, Hueschen and Richards (1379). Both faces of the 
germanium wafer were successively bombarded with B* ions at doses of 
10 14 lons/cm 2 at 25 keV and 2x20 14 lons/cm 2 at 50 keV. The wafers were 
next annealed at 20Q C C in an argon atmosphere to remove implantation 
damage and to activate the boron. The faces were then sputtered with 
500 A of palladium to serve as a wetting agent and finally a 1 urn layer 
of gold was sputtered on the faces for electrical contact. 

Individual detectors are sliced in 1 mm square sections from the 
wafer with a diamond wafer saw. The final crystals are 1 mm square by 3 
mm long along the [100] axis. The crystals are then etched in a solu- 
tion of HNO3, CH3CGOH, HF and HBr in the proportion 30: 18:18:1 to remove 
damage done by the diamond saw and any oxidized surface layers. During 
the etching process, the gold detector contacts must be protected with 
an overcoat of nail polish which is later removed with acetone. " 

Electrical contacts to the stressed detector are made through .002" 
thick soft copper pads which have fine insulated copper wires indium 
soldered to one corner. These pads are electrically Isolated from the 
detector housings by two layers of cigarette paper. The papers and 
copper pads are centered, one each, on the stainless steel plugs and 
held in place with fingernail polish* After the polish sets, excess 
paper may be scraped away with a razor blade leaving the Integrating 
cavity free of any non -metallic substances (save the detector). 


To construct a detector, the stainless steel plugs are first 
Inserted Into the Integrating cavity via the hole threaded for the set 
screw. The lead wires are pulled through to come out the entrance aper- 
ture. Next, the germanium chip Is Inserted into the housing through the 
entrance aperture and centered with a Teflon-restraining jig. Stress Is 
then applied through the set screw - ball bearing assembly. Application 
of enough stress to lower the detector resistance by 10-15% at room 
temperature is adequate to produce significant response out to - 195 
um. After application of the stress, the Teflon jig may be removed and 
the light cone mounted. The lead wires are carefully threaded out of 
the cavity through two tiny (#80 drill) holes. The detector will need 
to be restressed within a few hours of the initial stressing as the 
copper pads will cold flow somewhat. After this first restressing, 
however,, the detector will maintain a high level of stress throughout 
its useful life. (For some as yet unexplained reason, our stressed 
detectors seem to lose overall response on time scales of the order 6 
months « 1 year.) An advantage of these housings is that once detector 
crystals are cut and etched, a detector may be manufactured within an 
hour with a minimum of tools. This is an important point when one con- 
siders the rigors of shipping our spectrometer across the continent to 
NASA Ames Research Center. In fact, detectors have been made In the 
field and flown within a week of manufacture — though this was an emer- 
gency procedure not generally recommended. The manufacture of stressed 
detectors has been previously described by Watson (1982) and Haller, 
Hueschen and Richards (1979). 

For a more complete discussion on photoconductive devices see Bratt 




Data presented In this thesis were taken from two of NASA's air- 
borne observatories, the Lear Jet facility and the Kulper Airborne 
Observatory (KAO). 


2.5.1- Lear Jet. 

The Lear Jet Is by far the smaller of the two facilities. It 
carries’ a 30 cm f/6.5 Dali Klrklam Cassegrain telescope mounted In a 
side port. Due to space considerations, we are not able to mount 
our interferometer on the Lear Jet telescope and instead fly with just 
our grating spectrometer. The side looking dewar Is mounted directly to 
the backplate of the telescope which is Isolated from aircraft vibra- 
tions by an air-bearing and ayroscopically stabilized In flight to an 
RMS stability of - .5* on smooth flights. The air bearing requires 
that the cabin pressure be only - 5 Tb/in 51 higher than ambient pres- 
sure. Hence, the cabin altitude is maintained at - 14,000 feet and 
the experimenters and pilots must wear oxygen masks. The telescope 
looks out of a side port of the airplane at a relative bearing of 270° 
to the aircraft motion. The air bearing allows - 6° play In the tele- 
scope pointing In yaw and roll and the telescope may be mechanically 
cranked to center its window from 14° - 28° in azimuth. 
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There Is room for one experimenter on the aircraft. Unfortunately, 
data acquisition requires two experimenters. One experimenter, the 
tracker, sits in the back near the telescope. He operates the tracker 
rack and performs such tasks as balancing the telescope and cranking the 
elevation setting screw* He also has the heart-rending task of acquir- 
ing the field to be observed with the aid of a bore-sighted guide tele- 
scope. This telescope Is linked to a video display on the tracker rack 
through a photomultiplier tube and fiber optics bundle. The zoomed-in 
tracker field is roughly 1° in azimuth by .6° in elevation, and the 
faintest discernible stars are generally - 9 ^h magnitude on a 
moonless night. The zoomed out field is approximately ten times larger, 
and consequently it is difficult to see stars fainter than 4 ^h magni- 
tude zoomed out. The Lear Jet has automatic tracking ability facilita- 
ted by an on-board tracking computer. This computer senses the image of 
a guide star on the tracking screen and stabilizes that image through 
feedback to torqusr motors on the gyroscopes. If there are no guide 
stars bright enough to be tracked automatically (> 6^ or 7^ magni- 
tude), the observer may manually track with a joystick down to the 

dimmest visible stars. The chopper throw of the telescope (produced by 

* 

the oscillating secondary) is variable up to 14'. We typically use as 
large a throw as possible as our beam sizes are always greater than 

A A A 

4‘x6'; 6‘ along the direction of the chopper throw. We generally chop 
at a frequency which gives the best waveform and noise measurements: 
typically, We chop at 35 Hz. Chopper throw and beam size are determined 
through observations of point sources such as the planets. 
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The second observer spends his time running the experimenter's 
rack, primarily running the data acquisition program on our IHSAI 8080 
microcomputer* The data acquisition program was written by Scott D# 

Sny ers and recently revised and updated by Charles E. Fuller at Cornell 
Through Its use, the IMSAI sends pulses to the drive box which bursts 
the stepper motor controlling the grating and samples the signals Input 
to the A-D board. The oscillating signals from the spectrometer are 
synchronously demodulated by a pair of Ithaco lock-in amplifiers. The 
analog output from the Ithaco* s also goes to a strip chart on which the 
rack operator notes any unusual events (loss of track, high turbulence, 
etc.) In the data acquisition process. 

The optimal data acquisition technique is to make a measurement of 
a point In the plus beam, then in the minus beam, burst the grating and 
repeat the process until the entire number of points has been sampled in 
both beams. This mode, called plus-minus mode, would minimize the 
effects of temporal fluctuations in sky emissivity. However, nodding 
the telescope from plus beam to minus beam takes time and frequently 
results in loss of track on the Lear det. In addition, each single beam 
spectrum generally takes only - 90 seconds, so we usually decide to 
take single beam spectra, alternating minus and plus beams in successive 
spectra. That is, first a spectral run is taken with the source in the 
plus beam, then the grating is burst back to the starting position and a 
spectral run is taken with the source in the minus beam. Spectra gene- 
rally consist of - 12 points sampled, each point separated from the 
next by ~ 1/3 resolution element. In this manner, we sample sufficient 
continuum points about the expected line position to obtain good 
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calibration, and oversample the line position. The final spectrum, as 
discussed in Section 2.1 is the difference between the plus and minus 
beam spectral runs. Each run is recorded onto a floppy disk to be 
reduced post-flight. 

On the best flights up to 39 single beam spectra nave been taken 
even while mapping a region. Lear Jet flights are limited to roughly 
2^10^ in duration of which 65-70 minutes are available for astro- 
nomical observations. Whenever possible, we take a short (-10 min) 
leg on a calibrator and the rest of the available time on a source. 

This both verifies the position of the beams relative to the tracker 
screen and ensures similar atmospheric conditions for spectra taken on 
both objects. 

The Lear Jet, though a much smaller facility than the KAO offers 
several advantages that the KAO does not. For extended sources, the 
smaller telescope size of the Lear Jet does not matter as telescope 
aperture - solid angle product An (telescope aperture area x beam size), 
is a constant. Furthermore, the Lear chopper throw is significantly 
larger than the KAO chopper throw, thus it is much easier to find a 
blank sky reference when observing extended sources from the Lear Jet. 

The Lear telescope may observe closer to the horizon (14° vs, 35°) 
than the KAO* Thus, a larger portion of the sky is available to the 
Lear observer. One may observe the galactic center from the Lear from 
the latitude of Ames with no special flight planning. Lear flight 
series generally consist of 12-16 short (2^10' 11 ) flights while KAO 
series are generally 2-4 long (7^30 m ) flights. Thus, a Lear series 
allows one much more flexibility in flight planning — an important 
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consideration when timing Is Important as In occultation events (see 
Chapter 4). 

Finally, the Lear Jet typically flies at a somewhat higher altitude 
than the KAO. Most of the data taking leg of a Lear flight is flown at 
altitudes between 43,000 and 45,000 feet. While It Is true that the KAO 
may reach 45,000 feet, this flight level may be achieved only after ~ 

5 1/2 hours of the flight and 41,000 feet may not be achieved until * 

2 hrs into a flight. Thus, the Lear Jet will have advantages If the 
troposphere is high. 

2.5.2. KAO. 

The KAO is a Lockheed C141A Starlifter which has been modified to 
carry a 91.5 cm bent Cassegrain telescope. The telescope is isolated 
from the experimenter area in a cavity located just forward of the wing 
on the port side of the fuselage. As on the Lear Jet, the telescope is 
floated on an air bearing and gyrostabilized. However*, as the telescope 
is isolated from the experimenters' area by an ambient pressure cavity, 
the experimenters may work in a shirt sleeve environment with no oxygen 
masks. The cavity and telescope are usually pre-cool ed prior to takeoff 

»v 

to lower the thermal background. 

The KAO contains three on-board computers to control telescope 
pointing, data acquisition and data processing. NASA flight personnel 
typically consists of an ADAMS (Airborne Data Acquisition and Management 
System) operator who oversees KAO computer operations, a tracker opera- 
tor, who acquires the star field, one or two telescope operators and the 


flight director whose primary task Is to keep communication channels 
open between experimenters and all other personnel, especially pilots. 
There Is ample space on board for upwards of six experimenters, though 
we have generally flown with four or five. Generally, one person runs 
the data acquisition program, one writes notes on the strip chart and 
updates the rotation angle, a third reduces data on the EXEC computer, 
and the fourth "guides" which is to say tweaks the telescope beam to 
keep it centered in the spectrometer entrance aperture. 

The data acquisition program written by Charles Fuller runs on our 
IMSAI microcomputer which is linked to the ADAMS system through a serial 
port. The IMSAI controls the pertinent data acquisition commands: 
bursting the grating and interferometer, sampling the detector signals 
which have been synchronously demodulated by Ithaco lock-in amplifiers, 
commanding the tracker computer to nod the telescope through the DATA 
computer and receiving error flags from the EXEC through the DATA compu- 
ter. Data sampled by the IMSAI is written on a floppy disk and trans- 
ferred over to the ADAMS system where it is written on hard disk and 
magnetic tape. Data reduction is possible in flight using the program 
M4YAY, written by Scott Snyers. The program runs on the EXEC computer 
and performs fast Fourier transforms. 

KAO flights generally last 7 l/,2 hours, of which - 6 1/2 hours 
are available for observations. This is sufficient time to allow obser- 
vations of - six objects including a calibrator. Objects to be obser- 
ved must be within the 35°-72° window (in elevation) of the telescope. 

To ensure proper function of the interferometer, we usually observe a 
source with a strong continuum such as M42 or a planet on our first 
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leg. With a strong continuum source we single step our interferometer 
through zero path at two or three different grating positions to verify 
our preflight zero path position. Actual interferograms generally con- 
si st of 2 n +2 data points (n « 4,5,6) beginning with a sample one 
burst negative of zero path. For our 5 cm travel stage we have general- 
ly used 34 data points which fill a free spectral range large enough to 
encompass the bandpasses of both detectors. On the KAO,. nodding is 
usually a trivial procedure controlled by the IMSAI through the tracker 
computer. The nodding procedure employed generally consists of taking 
the integration at the first interferometer position in the plus beam, 
nodding and integrating in the minus beam, burst the Interferometer to 
position two, take a minus integration, nod to take a plus integration, 
burst and so on until we have completed the 34 integrations. Each data 
point typically consists of the difference between five seconds integra- 
ting in the plus beam and five seconds integrating in the minus beam. 
Including time spent in nodding (* 1 sec), pausing for telescope and 
electronic stabilization {~ 1 sec) and bursting the interferometer 
(" 1 sec) each 34 point interferogram takes - 8 minutes. Upon com- 
pletion of a forward interferogram, we take a reverse interferogram, 
beginning at largest travel and bursting back to zero path* During an 
entire flight over forty interferograms are often taken* 

Tracking on a guide star is accomplished by a tracking computer 
operated by NASA personnel. The tracker is similar to that on the Lear 
Jet, but much more sensitive, being able to track stars of visual magni- 
tude, my - 10*5 to a few arcseconds stability. Use of a dichroic to 
split the infrared beam from the visible enables us to guide on the 
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focal plane Image. We have made use of an excellent offset gulden 
developed by Prof, dames R. Houck to guide on stars as dim as 12^ 
magnitude in the focal plane. The offset guider is a twin reticle 
device which corrects for field rotation by actual rotation of one of 
the reticles. This movable reticle is placed over a guide star while 
the fixed reticle denotes the location of the source. The fixed reticle 
is aligned with the telescope ~ spectrometer boresight through a pre- 
flight procedure with a chopped~co1.11 mated infrared source. (Referred 
to as Erickson's source after its inventor* Dr. Edwin F» Erickson of 
NASA Ames Research Center.) Offsets to the guide star are calculated 
preflight and dialed into the movable reticle for each object, keeping 
the IR source in the entrance aperture is then just a matter of keeping 
the guide star in the rotating movable reticle. Rotation angles are 
provided from the housekeeping utilities of the ADAMS system. Focal 
plane guiding may be practiced on sources as far as 9‘ from their guide 
stars. Pointing accuracy may be quite high with such a system, as 

a 

small as ± 5", while guiding stability is probably of the same order 
RMS. The plate scale at the focal plane Is - 147mm which means the 

A A 

4x4 mm slits in our spectrometer define a square beam roughly l'xl 1 . 

A 

The chopper throw of the KAO is variable up to 9' , Generally we 
use a throw just larger than the source extent as nodding takes longer 
with larger chopper throws (it also takes longer with fainter guide 
stars) and the chopper is less efficient with larger throws (the time 
spent between beams is larger for larger throws). The KAO chopper is 
easily rotated up to 90°, and the centroid of the chopper duty cycle may 
be varied to minimize the infrared signal from the telescope. As in the 
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Lear Jet, chopper frequency Is chosen to minimize system noise. 
Generally chopper frequencies of 25 Hz have been used* 

2*6. Data Acqujs1ti<m-and Jtedudrt on* .. 

Data sampling rates and plus-minus reduction are essentially 
equivalent for the two observatories. We usually Integrate for five 
seconds in each beam for each interferometer or grating position. 

These five second Integrations, called points, P-j, are subdivided Into 
35 separate Integrations, s n , .143 seconds each, each of which is 
stored on floppy disk. A postflight da-glitching program written by C» 
Fuller, may be used to remove spikes in the data presumably caused by 
cosmic rays. They are easily recognized as spikes as they typically 
last one or two Ithaco time constants (usually set at 125 ms or - 1 Hz 
bandwidth). Generally data is automatically de-glitched by throwing out 
the five Integrations which deviate the most from the mean* This usual- 
ly results in " 20% increase in the signal /noise ratio. The data 
point for each beam is then given by the average of the N re- 
maining data points 

Pi v " VK iU 
n-l 

with noise, (standard deviation from the mean), given by 
m * t l/N I C (Pi - S„) 2 /(N-1) j i/2 

n=l .. 


63 


The fine!' sfgnal for each grating or Interferometer position Is the 
difference between tfte plus and minus beam measurements, so 

e - Pj w ; n\ « ( n 2 i* + n 2 -}* }^ 2 

thfcipint forward due to the Fourier transform routine required 
fcf* the KC data.,; .data reduction takes different routes for the two 
facilities. 

t>*r Jeta., 


?ratfcg spectra are averaged with the program AVRG. AVRG presents 
the jr:? with three options for data reduction. The common denominator 
fcr the three Is the final noise calculation. In all cases, two noise 
estimates are made, and the larger of the two Is chosen as the noise in 
the averaged file. If we average N files, each of which contains a 
measurement Pj of the source flux at wavelength A-j, then we define 
two error estimators 

1) Relative mean error — also known as the standard deviation from 
the meani 

«ras = ( 1/N Z C (P s - Pi) 2 /(M) ] ) 1/2 
1=1 

2) Mean error — often termed "internal error"} 

Oms “ ( 1/N l nj* ) lh 

If the number of spectra is sufficiently high, ( > 5 ) then the two 
error estimates should converge. 
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These two error estimates apply to a "straight 1 * average of the 
data* 1*e* where each data point Is given equal weight (« 1): 


h * i/« ! pi 

i«i 


If one feels that the errors associated with each data point Pj 
are reasonable, It Is desirable to perform a "weighted" average of the 
data files, where each data point is weighted by the Inverse of the 
variance: 


a j 1 P 1 /n 1 2 *j 1 1/ll i 2 

In this case, the mean error is given by; 

<*m * ( X Mij 2 ) 1/2 

The relative mean error is: 

*rmw “ l VN(N-l) l (fvP^/n-i 2 * J 1/nj 2 ] l/2 

If all of the n-j are equal, justifying the use of a straight average, 
then P w “ P s , c r nw * cfos and <? rfnw * cr rrns . 

External errors due to guiding, sky fluctuations or instrumental 
fluctuations may be compensated for through use of the third option of 
AVRGj maximum likelihood averaging of files. This technique was sugges- 
ted by Stephen E, Schneider and has been implemented for much of the 
Lear data reduction. The essence of the technique is that it acknow- 
ledges the existence of external systematic noise in our data, and tries 


to maximize the probability that the final average Is the correct one 
through adding a common external error, In quadrature to each data 
file. 

i 

Assuming the N data points Pf, for each wavelength are 
normally distributed with variance n^ 2 + s 2 where n-j is our measured 
noise and s is the external noise, to be determined, then the likelihood 
function is given bys 

L(s) e l//(2w) ,f 8 l/(ni 2 + s 2 ) 1 / 2 * 

N „ 

exp [-1/2 l ei 2 /(n^ 2 + s 2 )], 

1»l 

where ej 2 « P, and F « P w , as above. Maximizing with respect 

to s 2 we have 

(dL/Ss 2 ) | « 0 * 4/2 l l/(n-| 2 +s 2 ) 4/2 J 2 /(n-j 2 +s 2 ) 2 

^max * V 

If we let 

f(s 2 ) » l e| 2 /(n^ 2 +s 2 ) 2 « l l/(n-i 2 +s 2 ) « 0 
1 1 

then Sf{s 2 )/0s 2 * «2 l e-f 2 /(n-j 2 +s 2 ) 3 + l X/Cn^+s 2 ) 2 

i ■. i 

We may use the Newton-Raphson method to solve for s 

($‘) a » s 2 « f(s 2 )45f(s 2 )/ds 2 ] 

and iterate until | s 1 - sj $ an acceptable error. 
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The maximum likelihood option of AVRG makes these calculations to 
find s for each wavelength Aj. The new error for each file Is then 
nj ' ■ (n-j 2 * s 2 ) 1 ^ 2 and averages are calculated as In the weighted 
average options. The resulting average spectra differ little from 
weighted averages? means shift a few percent and the final errors are 
usually slightly smaller, as the relative mean error approaches the mean 
error. 

As mentioned In Section 2.1, spectra are divided by spectra of 
standard calibrators to remove as much as possible the convolved pre- 
sence of telluric water vapor features, detector and filter response 
functions, and beam profile from the source spectrum. Standard calibra- 
tors used on Lear Jet flights are the planets Venus, Mars, Saturn and 
Jupiter, and the Moon. Of the planets, Jupiter and Venus provide the 
most power, hence the best signal /noise ratios. However, the atmosphere 
of Jupiter is known to contain deep NH 3 features, and one must constant- 
ly be aware of such features when using Jupiter as a calibrator {see 
Figure 2.9). The J =* 2 and 0 a 3 rotation-inversion bands of NH 3 occur 
In the bandpass of our spectrometer at - 125 and 170 um (see Encrenaz 
et al. 1971 for details). Venus would be a more suitable candidate, but 
proximity to the sun prohibits Its use much of the time. The Moon might 
appear to be an ideal calibrator, however Its larger angular size means 
It fills the beam quite differently from the average astronomical 
source, which diminishes its value as a calibrator. 

Calibrator spectra, C-j, are taken at the same grating position, 
hence wavelength, Af, as the source spectra, Pj. The resulting 
spectrum s-j is then 


Sf “ P-j /C i 



Figure 2*9# Theoretical far-infrared spectrum of Jupiter from Encrenaz 
et al. 1971. The J - 2 (-170ym) and J * 3 (-125pm) pure rotational 
bands of NH 3 are within the spectral range of our spectrometer and 
must be avoided when using Jupiter as a calibrator. 
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with associated noise cr^ * $•{ (acij/Ci 2 + n l 2 / p l a ) 1 ^ 2 * Absolute flux 
calibration follows from the known physical parameters of the cali- 
brator. 

Dividing of spectra is accomplished by the program UNFIT. UNFIT 
also allows one to calculate the flux in an emission or absorption 
feature by fitting a weighted least squares Gaussian plus polynomial to 
the divided spectrum. 8asel1ne curvature is usually entirely removed 
through the dividing of spectra, so the polynomial is usually first 
order. 

If one knows the precise location of an expected line in low 
signal/noise ratio data, one may recover some information from the data 
through the binning option of UNFIT. Since we oversample our data, 
data, usually three samples/resolution element, we may average, or bln, 
these three elements into one number. This final number generally has a 
S/H ratio /3 times better than tho individual samples, resulting In the 
recovery of useful data front a noisy spectrum. Binning may be achieved 
through a straight or weighted average, as above. 


2.6,2. KAO. 

The data acquisition program written by Charles E# Fuller performs 
in-flight plus -minus reduction on the KAO, This plus-minus reduced file 
is transferred to the EXEC via the DATA computer. The data file is 
accessible to the program, M4YAY, written by S. Snyers which allows the 
experimenter to perform a first cut data reduction In flight. The fast 
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Fourier transform routine for M4YAY was developed by B. Langdon and S. 
Sande at Princeton University, from the approach of Cooley and Tukey 
(1965) . In-flight data reduction Is an Invaluable aid, as it enables 
us to estimate when we have Integrated sufficiently on a line and may 
move on to the next feature. 

The KAO data Is de-gl Itched prior to post-flight plus-minus reduc- 
tion, using the Individual Integrations stored on floppy disks during 
flight. The Fourier transform program which runs on the IMSAI, again 
written by S. Snyers, Is similar to the in-flight program, The program 
reconstructs wavelengths from the burst size, step size and grating 
central wavelength as In Section 2,3.3 above B The Fourier transform 
routine calculates both the cosine and sine transforms, As reflected 
data are symmetric about zero path, all the sine components are zero. 
Averages are performed in the standard way, with each Interferogram 
being transformed first and resultant spectra averaged with equal 
weight, i»e. 

Sj * 1/N l Sij 5 nj * [ 1/N l Csij"Sj)V{N~l) j l/2 
1 =1 I s ! 

where Sj « average of the N data points s-jj each with error n-j j . 

The noise nj is just the standard deviation from the mean. 

With our 5 cm travel stage, walk off has been a chronic problem. As 
mentioned in section 2.3,5 above, this "droopy shoulder" appearance may 
be minimized by fitting a line to the data and raising the wings. In 
addition, our Fourier transform programs allow one to apodize the data 


with spectral lag windows. These windows smoothly apodlze the data from 
a weight of one at zero path to a weight of zero at largest path. 'This 
effectively eliminates the discontinuity at largest path and results In 
a higher signal to noise ratio In the transformed spectrum at the 
expense of system resolution. When apodlzatlon Is necessary we have 
generally used the Bartlett (triangle) window which weights the data 
according to: 

w i * 1 * 1/N 

where 1 denotes the data point as measured from the zeroth data point 
(zero path) and N is the total number of bursts (* # data points - 1). 
The CO spectra of Figures 2.1 and 2.4 are apodized data with droopy 
shoulders removed. For a complete discussion of spectral lag windows, 
see Jenkins and Watts (1968), For a complete discussion of the fast 
Fourier transform see Bergland (1969), A discussion of Fourier trans- 
form spectrometers and 1 nt erf ero gram sampling error correction tech* 
niques may be found In Holler and Rothshild (1971), 


CHAPTER 3 


FINE STRUCTURE LINES AND RADIATIVE TRANSFER 


3 .1._ Introductions 


This chapter Is split Into two sections. The first section dls « 
cusses the physics of fine structure transitions and radiative transfer, 
while the second section briefly discusses some of the processes which 
compete with fine structure lines to cool the Interstellar medium. 

Fine structure levels occur In any atomic or Ionic species with 
valence electrons. The two species of concern in this work, CII and 01, 
have one and four equivalent 2p electrons respectively. The resulting 
2p x and 2^ configurations are ordered by the residual electrostatic 
interaction into terms which differ in total orbital angular momentum 
L = X If and total spin angular momentum of the electrons S s l Sf. 

We have denoted by If the orbital angular momentum, (| If j “ 1), 
and by s-j the spin angular momentum (| Sf j - 1/2) of each 2p elec** 
tron. For CII there is only one 2p electron resulting in the single 
term 2 P. The superscript denotes the number of levels of which the term 
is composed, and equals 2S+1. The letter P denotes the total orbital 
angular momentum of the term and follows the usual convention S,P,D,F, 
... refering to terms of L - 0,1, 2, 3,.., , The 2 P term is further split 


by the spin-orbit (L*S coupling) Interaction Into the levels 2 p 3 / 2 > 

2p i /2 where the subscript denotes Os L + $* Each level* d, has a 
degeneracy gj given by gj * 2dH. The 157 pm CCIID fine struc- 
ture line results from radiative transitions from the 2 P 3/2 level to 
the 2 Pi /2 ground state# 

The configuration of 01 is split into the three terms l S, X D, 
and 3 P. The lowest lying term is that with the highest spin (Hund's 
rule), the 3 P term. This triplet is split into the levels 3 Pq» 3p i and 
3 P 2 * As the Sp shell is more than half full, the levels are inverted, 
with the level of largest d, 3 P 2 , being the ground state. Figure 3.1 
shows the 01 level diagram. The l S and l D levels are metastable, the 
forbidden transitions to the 3 P term resulting in the familiar nebular 
lines at - 6300 A and ~ 2972 A. These levels have excitation tern-" 
peratures 2 20,000°K and are thus negligibly excited in the low tempera- 
ture (T ~ a few hundred °K) regions where the FIR fine structure lines 
dominate the cooling. The 145.53 pm Tine arises from radiative transi- 
tions from the 3 Pq - 3 Pi level while the 63.17 pm line occurs through 
transitions from the 3 Pi level to the 3 P 2 ground state. 

Fine structure lines occur between levels in a configuration. The 
fine structure lines of concern in this work occur between levels within 
a term and as such involve no change in the spatial wavefunction of the 
electronic state, but rather changes in d, i.e», the vector addition of 
the S and L vectors. As such, these transitions are forbidden to 
electric dipole radiation. Magnetic dipole transitions are allowed and 
hence the excited fine structure levels are metastable with lifetimes 


48,60Q°K 



Figure 3.1. Energy levels for the (ls) 2 (2s) 2 (2p) 4 configuration of 
neutral oxygen. Wavelengths for the visible and ultraviolet lines are 
taken from Osterbrock, 1974. 
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Much may be learned from the study of FIR fine structure lines. As 
we show In section 3.3 below, fine structure lines, particularly the 157 
urn Cell] line dominate the cooling of gas at temperatures below a few 
thousand degrees. In section 3.2, we show that the measurement of one 
FIR fine structure line together with a priori knowledge of the number 
density of exciting particles, n c {such as from radio measurements) 
allows one to determine the relative abundance of the emitting species. 
The ratio of the two CO!] fine structure line Intensities enables one to 
determine n c In a model -Independent way. The Importance of far 
infrared and submlllimater fine structure lines has been the subject of 
several reviews in the recent literature including: Simpson (1975), 
Kafatos and Lynch (1980), Watson and Storey (1980), and Harwlt (1984). 

M & att ja Jcan g f 

3.2.1. . Equation of Transfer. 

The fundamental equation governing macroscopic radiative processes 
in regions as diverse as stellar interiors and diffuse interstellar 
clouds is the equation of transfer: 

' dly/dS * -Ctyly + j v [3,1] 

where I v [erg s wl cm** 2 sr" 1 Hz” 1 ] is the specific intensity of the 
beam along a pathlength s [cm], and ct v [cm“ l ] and j v [erg s" 1 sr" 1 
cm“ 3 Hz” 1 ] are the monochromatic absorption and emission coefficients. 
This equation fully describes the macroscopic behavior of radiation in a 
medium which both absorbs and emits radiation. The measured quantity, 
Iy, allows us to determine the physical properties ay and j v of 


the medium. The primary task for this section Is to determine forms for 
these two coefficients* and hence solve 34. Much of the formalism we 

t. 

use in this discussion follows Ryblcki and Llghtman (1979) and Spltzer 
(1978). 

In principle, when c v and j v are known it Is easy to solve the 
equation of transfer. Unfortunately, scattering is also present in a 
real medium thus complicating the analysis. We postpone the discussion 
of scattering to section 3,2.4. 

Equation 34 is readily solved In two limiting cases: 

1 Emission only, a * p. If MP) is the background intensity 
Incident upon an emitting slab of thickness s, then the emergent 
intensity I v is given by 

I v */ Ms') ds‘ + I V (P) 
o 

2 — Absorption only, j v * P» Let I V (P) again be the intensity 
incident upon our absorbing slab of thickness s, then the intensity will 
be attenuated In the medium so that the emergent Intensity is given by: 

I v “ M0) exp('**J ct(s') ds * ) 
o 

One defines the optical depth, t v , by 

s 

dT V “ cfy ds + t v (s) “ / Oyfs 1 ) ds 1 . [3.2] 

o 

Hence, in case 2, 

I. - MP) e“ Tv 
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From the exponential form of the above equation, It Is clear that 
the average photon will travel to (at least) an optical depth t v with 
probability e’’ Tv * Thus, the average optical depth traversed by a 
photon, «r v >, Is unity. 


w 


<t v > 5 / t v 6 hTv dr v 3 1 
0 

We may define the mean free path of a photon, Ay, In a homogeneous 
medium by <t v > 3 <a v Ay> B 1 Ay s l/a v . A medium Is said 
to be optically thin if the average photon escapes the medium without 
absorption, i.e., t v < 1* If t v > 1, the medium is optically 
thick. Dividing 3,1 by ct v we may rewrite the equation of transfer in 
terms of the source function, S v s j v /a v and optical depth t v : 


dl 

dr. 


v _ 


"Iy + Sy 


This leads us to the formal solution of the equation of transfer: 


ly(Ty) 3 Iy(P)e^V + /Y (t V" T v')Sy(T v ) dty 1 [3.3] 


If the source function is constant in the medium, we may integrate 
to yield 

Vfrv) * s v + 

Note: As t v + «, I v + S v , Also note that in thermodynamic 
equilibrium, we must have j v 3 a v B v (Kirchoff 's law) where B v is 
the Planck bl ackbody function given by 
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B, ■ ( e llv/kT -l )-» 

c 2 

Thus If a system Is thermal Ised at some temperature T, the source 
function reduces to the Planck function. 

3. 2^2* Einstein Coefficients and Detailed Balance. , 

Einstein first realized that Kfrchoff’s law Implies a relationship 
between the emission and absorption of radiation on the atomic level. 
This relationship is expresses through the Einstein coefficients* 

Let us begin by limiting our attention to a two level quantum 
system, denoted by u (upper) and 1 (lower). The difference in energy 
between the levels is aE u i B hv u -) • We define A U | as the sponta- 
neous decay rate through emission of a photon when the system makes a 
spontaneous transition from level u to 1. In addition, a passing photon 
of energy hv U; j may be absorbed by the system in the ground state 1 and 
excite the atom to the excited level u* If we denote 0 V * l/4ir ^I v 
dn as the mean intensity of the photon field, then the probability, 

P a , of a photon being absorbed per unit time is given by 

P a “ B J du 

0 

where B| U is the Einstein coefficient for absorption and <j>(v) is the 


11ns profile, which in general Is strongly peaked about v u } and Is 

C4 

normalized by / $(v) dv * 1. 

o 

If the system Is In the excited state, then the photon field may 
stimulate emission of a photon* The Einstein coefficient for stimulated 
emission B u j Is defined such that the probability, P s , per unit time 
of the emission of a photon through stimulated emission is.* 

p s s B ul/ tyW dv 
0 

For Individual atoms, the line profiles may be approximated by a 
Lorentz profile: 



Y/4tt 2 

(v«v u -|) 2 + (y/4tt) z 


where y is the natural line width, related by the uncertainty principle 
to the Einstein A coefficient for the transition by y ** A u :j. For 
the lines of interest In this work, A u i * lb" 1 * -IP** 6 s wl . At reason** 
able astrophysical temperatures (T > 10°K), Doppler line widths are 
M’Dop > 10 s «10 7 Hz for these lines — much greater than the natural line 
widths. Therefore, the mean Intensity J v is constant over the natural 
line profile, and we may take 0 V out of the integral above. Eval.ua- 
ting, at the line center, we have: 





C3.43 


The Einstein coefficients are related at the atomic level (the 
system need not be in thermodynamic equilibrium) by: 
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B lu s V9l B ul 

£3.53 

c 2 

£3.63 

ul 2hv^ ul 


The absorption and emission coefficients may thus be written In 
terms of the Einstein coefficients: 

« v « hv/4ir ( n-jB^ - n u B ul ) <{>(v) [3.7] 

J v « hv/4ir ( n u A ul ) <f.(v) [3.83 

Stimulated emission has been included as negative absorption since stim- 
ulated emission and absorption both are dependent on the local photon 
field. 

Rewriting the equation of transfer 3.1 in terms of the Einstein 
coefficients we have: 

~ ( n 1 B lu - n u E ul H(v)I v + ^ n„A ul4 ,(v) C3.9] 

We divide 3.9 by a v and integrate to obtain: 

*v " o n ^l e " Vds [3.103 

■.» 

Most far-infrared fine structure lines are optically thin (t v < 1). 

For these lines, we may immediately integrate 3,10 to yield: 

« hv hv n u 

1=1 I u dv = — N u A ul = Hj A u1 £3413 

0 Air 4 it rij 
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Thus our measured parameter, the intensity in the line Integrated over 
the tine profile, I, determines the column density N u s n u s of the 
emitting species In the excited state u* To determine the total column 
density, Nj of the emitting species, we must determine the relative 
level populations n u /ni* We turn our attention to an atomic system 
with arbitrary number of energy levels nw 

In local thermodynamic equilibrium (LTE), levels are populated 
according to the Boltzmann distribution: 

ly/nj * gj/f e" E i /kT [3.12] 

where gj is the statistical weight of level j which has an excitation 
energy Ej above the ground state* nj/nj is the ratio of the 
number density of level j to the total number density of the species 
nj. f is the partition function defined by: 

f »? 9 k e-VkT 
k=9 

. . * 

LTE is maintained only when collisions! processes dominate radiative 
processes in determining the level populations* In the interstellar 
medium, where number densities are quite small, colllsional excitation 
rates for fine structure levels are often comparable to radiative tran~ 
sition rates* Thus, the fine structure levels are rarely populated in 
accordance with the Boltzmann distribution. For a steady state system, 
we may determine the* relative level populations through the principle of 
detailed balance: the number of transitions into a level j per unit 
time must equal the number of transitions out of that level per unit 
time* 
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Levels are populated and depopulated through both collislonal and 
radiative processes* For a steady state, the net change In the popula- 
tion of a level j per unit time Is zero, dnj/dt * 0* Equating transi- 
tions into a level 4 to those out of the level j we haves 

Ik (leWijt Vv w ) ^ %% 

* n jU>cTJk'%‘W + Z .Ajk3 C3-133 

Equation 3.13 Is subject to the constraint that l nj * nj — 
the sum of the populations of the levels equals the total population of 
the emitting species. In 3.13 above, n k Ay (njAj k ) are the 
spontaneous decay rates Into (out of) the leva! j, and 
(njBj}jd V j k ) are the absorption and stimulated emission rates 
Into (out of) the level j, d V jg Is the mean Intensity of the 
radiation field at the line center frequency v k j. The colli si onal 
excitation rate Into the love! j Is defined by n c v k j Csec Ml 3 where 
again n c « number density of collislonally exciting particles 
(electrons, protons, hydrogen atoms or hydrogen molecules), and Y k j 
Is the rate coefficient for collislonal excitation of level j. The rate 
coefficient for collislonal de-population of level j, Yj k is related 
to if kj bys 

Vjk = Vj *kj e “ Ej ' kAT 

where Ej k = E k - Ej. 

Rate coefficients are defined as the mean value of the transition 
probability per unit time for excitation of a level averaged over a 
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Maxwellian velocity distribution for exciting particles. The colllsonal 
excitation rate Is given by * <v c ojk> where <7j|> Is the 
excitation cross section, v c Is the velocity of the exciting particles 
and the brackets denote the velocity average. For electronic excitation 
of ionic species* the cross section may be written in a particularly 
useful form by defining a "coll Ison strength'* ft(j,k): 

ffjk B w /9j Ch/(2ffjn Q v e )3 2 n(j,k) 

where % « electron mass. Collision strengths are calculated quantum 
mechanically and are fairly constant with temperature above the 
threshold energy for colllsional excitation of the ion. For the 
fi ne-structure transitions of Interest In this thesis, the electron 
temperature In Ionized regions is always much greater than the 
excitation threshold. Collision strengths may be shown to be symmetric, 
n(k,j) « n(j,k) and are related to the rate coefficients by: 

h 2 n(j,k) 

^ 9 k ( 27?!!^ ) 3 7 2 ( kT ) 1 7 2 ' 

Thus above threshold, the rate coefficients for electronic excitation of 
ionic species are proportional to T wl/ ^ 2 . 

Colllsional cross sections, <?kj, tend to be constant for low 
energy collisions between neutral particles, such as HI - 01 
collisions. Collision strengths thus vary as T resulting in rate 
coefficients proportional to T 4 * 1 ^ 2 . Rate coefficients must thus be 
tabulated as a function of temperature. 
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Table 5,2 
Observations 


Region 

Transition 

Intensity 
(erg s" l cnr 2 sr“ l ) 

Comments 

Broadbeam 

(W) 

[01] 3 Po- 3 Pl 
(145,53 pm) 

2,240, 4xHT 3 

This work normalized 

A A 

from 7* x7* beam. 


[01] 3 p 1 - 3 p 2 

(63.17 nm) 

3.0*10" 2 

Estimate from contour map 
of Furniss et at. (1983), 


[Oil] 2 P 3/2 - 2 
(157.73um) 

Pi /2 3,84l.lxl0" 3 

Kurtz et al. (1983), 

Trapezium 

a a 

(1W) 

CO 13 3 Po~ 3 Pi 
(145.53 pm) 

5,741,7x 10 m3 

This work. 


[01] 3 Pi- 3 Pj 
(63.17 |int) 

6,240, 6xl0" 2 

Werner et al, (1982) 

A 

(45“ circular beam). 


Cr.3 


5»4.JL Ionization Fronts,. 


The 63 ynv Intensity expected from an ionization front of thickness 
D and number density of electrons, n e is given by equation 3,19; 

T j63 / 2 n\ n 0 

I - (ne* D) eqt — 

n e n 0I n H 

In this equation, and for the discussion which follows we use gaseous 
oxygen and carbon abundances as determined for the Inner regions (within 
V of the Trapezium) by Peimbert (1982); e 0 n 0 /n H - 4.7x10"^; 
e C n C^ n H * 3.7x10"**. 

The thickness of the ionization front is likely to be of the order 
of the mean free path of a Lyman a photon, D - 10 17 /n^ cm. If we 
assume pressure equilibrium between the D-type ionization front (cf. 
Spitzer 1978) and the inner regions of the Orion Nebula [n e - 4500 
cm" 3 ; T e - 8900°K (Peimbert 1982)], 50% ionization and an electron 
temperature T e - 6000°K in the ionization front, then; 

I * 4xl0" 6 erg s" 1 cm" 2 sr" 1 

where we have used Figure 3.7 for j6$/(ne n oi)* This is four 
orders of magnitude smaller than the observed [01] intensities. We 
conclude that the observed [01] emission does not arise in the ioniza- 
tion front. 


5 0 4»2. Warm Neutral Regions* 


C0I3 transitions may also be excited through collisions with neu- 
tral hydrogen atoms In the dense CII regions just beyond the Ionization 
front. (See Chapter 4 for a discussion of CII regions.) We begin our 
discussion by first restricting ourselves to our broad-beam (4‘x6' emit- 
ting region) measurements. 


5,4*2,!, Broadbeam models. 

CII regions typically have temperatures of the order SQO’K, In 
pressure equilibrium, the CII region would thus have a number density of 
HI atoms nuj ■ 3xl0 5 cm" 3 . For such temperatures and pressures, we 
would expect a line ratio, Ig 3 I 4 S ~ 40 (Fig* 3,8), However, the 
observed intensity ratio for our broadbeam measurements is - 18, This 
means the 63 urn radiation must be self-absorbed. 

To reproduce the observed 63 pm intensity, the minimum column den- 
sity of neutral oxygen in an emitting region of temperature T - 300°K 
and density n^ " 3xX0 5 cm" 3 is given by equation 3,19s 

No i = — * 1.6xl0 18 cm" 2 

. (as3/n H n 0I )n H 

The line center optical depth is given by eq. 3.22 in equation 3.24: 
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where N 2 and N* are the column densities of oxygen atoms in levels 0*2 
and 0*1, e 6 3 is the probability per absorption of colli sional 
de-excitation of the 0*1 level, gj are the statistical weights of 
level 0, and Av D is the Doppler width of the 63 m line. For T * 
300°K and n HI * 3xlU 5 cm" 3 , e 63 - .38 (equations 3,26 and 3,27). 

Using Figure 3,4 we estimate the relative level populations: nx/r^ - 
,052, no/n 2 - .019. Thus N 2 - .93 Hqj , - .052 Nqi, and we have: 


3xl0 5 



Oaffe and Pankonin (1978) have measured the C109a and C137& radio 
recombination lines from carbon atoms in the dense CII regions just 
beyond the ionization front. They measure a line width (FWHM) M 4 km 
s" 1 , corresponding to a Doppler width, Av^ (see equation 3,23), of 
2,4 km s" 1 . If we assume that this is the width of FIR lines emerging 
from the dense CII regions, the 63 ym line must have^ an optical depth 
*63 £ 1»2, We construct models for the [01] emitting CII regions with 
self-absorption of the 63 ym radiation in mind. 

Radial velocities of the carbon recombination lines led Jaffe and 
Pankonin to model the emission as arising from three distinct dense CII 
regions, distinguished by their radial velocities. The low velocity 
( v lsr " ® km s" 5, ) CII region lies in front of the HII region in the 
dark bay. The middle velocity region (vlsr ** 8 *5 km s" 1 ) is thought 
to be a sheet of material viewed edge-on which lies at the interface 
between the HII region and M43 to the northeast. The highest velocity 
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V LSR * ^ 5-1 ) Ctl region lies behind the Trapezium at the Inter- 

face between the HII region and the molecular cloud. The half -power 
distribution of these CII regions is * V (R.A.) x 5' (Dec.) in extent 
— centered roughly 2* north of the Trapezium. Thus the CII regions 
coincided reasonably well with the observed [01] emission. 

Kurtz et al. (1983) have observed the 157 urn [CII] line emission 
from these CII regions using an earlier version of our grating spectro- 

A A 

meter. In their 4 , x7 l beam centered on the FIR continuum peak, they 
measure a line intensity Iis 7 " 3,8xlCT 3 erg s*" 1 cm** 2 sr~ l # Assuming 
the 145 um [01] line and the 157 urn [CII] lines arise in the same 
regions, we may use the ratio of the line intensities (assuming both 
lines are optically thin) to determine the temperature of the emitting 
regions. Figure 5.3 displays the ratio of the reduced emissivities of 
the two lines: 

f ^ 145 10 i 

R 2 

( Jlg7 \ ) 

n CII n Hl n H 

as a function of density for various temperatures of interest. The 
cosmic abundance terms above (n 0 /n^ * 4.5xlp‘ ,tf ; n c /n^ * 3.7xlO M ' f ) 
normalize the emissivities to emissivity per hydrogen atom. 

Also plotted in Figure 5.3 Is the locus of points determined by the 
requirement that there be pressure equilibrium between the compact HII 
region and the CII fronts. It can be seen that the 145 pm to 157 pm 
line ratio is a sensitive temperature indicator in the density regions 
of interest. The ratio of the measured 145 pm intensity to the 157 pm 
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t 


* 


Intensity is: 

r = ii!il « q.58 

*157 

Using Figure 5,3 we determine T ** 350% n^ - 2,3x10 s cm" 3 for the 
emitting regions. 

From equation 3,19, we calculate the average line of sight column 
density of [013 emitting material*. 



*145 


{Jiffs/ r, oi n HI^ n 0^ n H ^ n HI^ 


« 7»3xl0 21 cm** 2 


where we have used Figure 3,6 to determine for T « 

350°K and n HI * 2.3xl0 5 ctrr 3 . This column density is within the 
range of models for CII regions (see section 4,3.2,1), 

Using equation 3,22 In equation 3,24, we calculate the line center 
optical depth of the 145 pm line: 


t 14S 


s}(|a 3 Aii*5 


Av 8rr 


372 


C22.ui 

a 



may be calculated from equations 3.28 In equation 3.29, evaluating 
0 V with equation 3,44: 

e l45 * *76 

Using Figure 3,4 to obtain the relative level populations we have: 

t 145 ~ *64 


In a completely analogous manner, we calculate e 157 * *98, n 3 / 2 /ni /2 

- 1*55, N 3/2 * 1.62xlG 18 cm" 2 , W\jz " 1.08xl0 18 cm" 2 , so that T 157 

- .34. Thus both lines are optically thin, and we are justified In 
using their ratio to determine the temperature of the emitting regions. 

As argued above, the 63 pm radiation must be self-absorbed. To 
estimate the 63 pm Intensity from the CII regions, we must refer to the 
analysis of optically thick lines given In Section 3.2.5 above. 

The Intensity of the 63 pm line as related to the Planck blackbody 
function, By Is given by equation 3,43: 

I S9 = IyAv = B V 4 V [ 1 * a - l 1 — ■ - ] [5.2] 

( e B„ e -S) (e&-e"&) 

where By Is the Planck blackbody formula, a a (l-/e)/(X+/e') and £ 4 
/3e t 63 , and av [Hz] Is the line width, full width half maximum. Having 
solved for T and n^j above, s 63 is determined. Using equations 3.27 
and 3.44 in 3.26, for I 33 * .030 erg s " 1 cm " 2 sr" l » we obtain s 33 s 
,38, The only free parameter in 5.2 above is xe 3 , the line center opti- 
cal depth, given by equation 5,1. Using Figure 3,4 to find the relative 
level populations at T « 350% n^j a 2.3x10 s cm " 3 we have: 

3.9 

*63 - 

AV D 

where Av d is expressed in km s~ l and we have used Nqi ~.4,5x10~ 4 
Nhi S ~ 7.3xl0 21 cm" 2 . To produce the observed 63 pm intensity 
(I 63 - ,030 erg s" 1 cm" 2 sr" 1 ) we require t 63 - 2,0. The Doppler width 
of the 1 i n e i s th en Av^ * X, 9 km s" 1 , or a line width (FWHM) of 
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Av * 3.2 km s“ l 

This line width Is somewhat smaller than, but consistent with the C109a 
recombination line widths measured by daffe and Pankonln. Their mea* 

sured line widths range from 1.7-7.6 km s'* 1 with an average line width 

* 

of 3,9tl,2 km s'* 1 , Thus, the C0I3 lines most probably arise from HI 
Impact excitation of 01 In the carbon recombination line emitting CII 
regl on. 


5.4v2.2* Trapezium region model# 

To model the enhanced emission from the Trapezium region, we must 
first subtract the contribution due to the broadbeam model. This leaves 
us with [013 intensities: 

I S3 « ,032 j I 145 « 3.5*1. 7 x 10 m3 erg s“ l cm” 2 sr" 1 

The optical sodium D line work of Isobe (1980) and the b1dimen« 
sional photometry of Uques and Vidal (1979) reveal the presence of 
globules, i.e. neutral gas condensations embedded within the Trapezium 
cluster. These globules could be remnants of the original molecular 
cloud from which the Trapezium stars formed [evaporation timescales - 
3x10** years, consistent with the age of the cluster (Vidal, 1982)] or 
might have formed from elephant trunk Intrusions of neutral gas into the 
HII region, a la Pottasch (1965), 

Each globule (typical diameters $ .002 pc) will have an associated 


CII front enveloping It, in much the same way CII regions envelop GHC's 
In the galactic plane (Chapter 4). In many cases the carbon will be 
Ionized throughout the globule then the entire neutral region becomes 
a good source of [01] emission* 

If we assume the conditions In these CII regions to be similar to 
those from which our broadbeam emission arises (T - 350°K, - 

2*3xlQ 5 cm” 3 ) and take a line width of “ 4 km s* 4 , then the measured 
63 pm Intensity requires a line center optical depth 733 - 2.0* This 
corresponds to a column density of atomic hydrogen in the [01] emitting 
regions of the order Nm ** 10 22 cm” 2 . The optically thin 145 pm 
intensity is then: 

I^ 5 a 3,0xl0** 3 erg s” 1 cm” 2 sr” 1 

This is smaller than the observed intensity but well within the 
measurement errors. 

We thus require roughly one globule along a typical line of sight 
in the Trapezium to explain the observed emission. This is consistent 
with the conclusions of Isobe (1980), The extra CII region in the 
Trapezium would also give rise to enhanced carbon recombination line 
emission, which is also measured. Our results are summarized in Table 


Table 5.3 
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We have made the first astronomical detections of the 145,53 um 
C0I3 line through observations of the Orion Nebula* By comparing 
observations of the 3 P* ** 3 Po (157,74 um) transition of Cl! we are able 
to deduce that these three emission lines all arise In dense CII regions 
at the Interfaces between the Ionized gas and molecular cloud In Orion, 
Using the 145,53 um to 157,73 urn line Intensity ratios, we are able to 
determine that these regions are at a temperature T * 350°K with 
number density of hydrogen atoms n^i * 2.3x10 s cm” 3 . We show that 
the 63,17 urn COI3 line is self-absorbed with line center optical depth 
t 63 * 2.0. This is the first such determination of an optically thick 
FIR fine structure Tine. 


APPENDIX 


Selected Far* Infrared Fine Structure Transitions % 20 i«i) 
(After Schntid-Sur^, 1982) 


Transition 

Species 

Excitation 

Potential 3 

(eV) 

Ionization 

Potential 3 

(eV) 

Wavelength 

(microns) 

First 

Astronomical 

Observation 

2 P3~ 2 Pl 

M *4 

BI 

HH 

8.298 

655,57 

HHH 

2 2 

CII 

11,250 

24.383 

167,737i.002b 

Russell et al. (1980) 


win 

29.601 

47.448 

57.343i-.003 

Moorwood et al. (1980b) 

f 

OIV 

54.934 

77.412 

25.87 ±.02 

Forrest, McCarthy, and 
Houck (1980) 


All 


5.986 

89.237^.0025 

NHN 


sin 

8,151 

16.345 

34,814 

MHH 


pin 

19.725 

30.18 

17,885+ .007 

HHH 

2 Pl« 2 p3 

H « 

2 2 

FI 

MM 

17,422 

24.745±,001 

MMM 




192 


Excitation 

Potentials 

Ionization 

Potential^ 

IHHSH 

MM 

11.260 

14,534 

29,601 

35,116 

54.934 

62,707 

87,138 

97.11 

126.21 

MH 

8,151 

10,486 

19,725 

23,33 

34.83 

39.61 

53,46 


51,816+. 003 
88.355+.008 

25,83 

44,25 

14,32 

24,28 ±,02 

68,473 

129.68 

32.87 

60,68 

18.7129 
i 0004 
33.456+. 001 

11.76 +.01 
20,38 ±.02 


Astronomical 

Observation 


Oaffe et al« (1984) 
Phillips et al. (1980) 


Melnick et al. (1978) 
Ward et al, (1975) 


Forrest, McCarthy, and 
Houck (1980) 


Baluteau et al. (1976) 
Moorwood et al. (1980a) 







TranslHlonl Sped es I Excitation | Ionlzati onl Wavelength . 

Potential 3 1 Potential I i Astronomical 

Observation 


“ 


FII 17.422 
Neill 40.962 
NalV 71.64 
SI 

CUI 12.967 
Aril I 27.629 


s D 3 « 5 D 4 Fel 

Fell! 16.18 

5 “2- S[ >3 Fel 

Felll 

5n. J3n- CoT 


Di« s D 2 Fel 


13,618 

145,52547 

£.000085 

63,17000 

£.00003 

34.970 

67.16 

29,33 

63,45 

36.02 

15.55 

98.91 

21,29 

9.039 

10.360 

56.309£»01 

25,2490 

£.0003 

23.81 

33,281 

14.368 

40,74 

21,84£,1 

8,9910 

£.0002 


7.87 

30.65 


Fel 

Felll 


24.04 
22.93 

34.71 

33.04 

54,31 

51,68 

111.18 

105.4 


Shure et al. (1984) 
Pottasch et al. (1984) 


SHure et al, (19S4a) c 
Gillett, Forrest, and 
Merrill (1973) 
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Transition 

Species 

Excitation 
Potential a 
CeV) 

Ionization 

Potential 3 

(eV) 

Wavelength 

(microns) 

First 

Astronomical 

Observation 

%»% 

Fell 

•7.870 

16*18 

25*99 

MHNi . 

2 2 










35*352 


2 2 






6 d 3 - 6 d 5 




51*28 


2 2 






6 DiH 6 D 3 




87*41 

Hhw 

2 2 






V 5 d 3 

FeV 

54.8 

75.0 

20,79 


5 d 3 - 5 d 2 




25,97 


5 o 2 - 5 Di 




36*50 

(MM*# 

- 

s o r 5 o 0 




68.97 



a From Lang (1980) * 

Excitation Potential s Energy needed to form the species. 
Ionization Potential 5 Energy needed to ionize the species, 

k Crawford et al , (1985) , 

c Possible detection. 


REFERENCES 


* 

Allen, C,W« (1975) Astrophys! cal Quantities, 3 rc * Ed* (Londons Athlone 
Press}* 

Altenhoff, W*J*j Downes, D., Goad, L., Maxwell, A*, and Rinehart, R* 
(1970) Astron . Astrophys , Suppl , Ser , 1* 319, 

Armstrong, K.R., and Low, F.O, (1973) Appl ♦ Opt , 13, 2007. 

Bahcall, J.N., and Wolf, R.A. (1968) Astrophys . J, 152, 701* 

Baluteau, J,~P,, Bussoletti, E,, Anderegg, M., Moorwood, A.F.M., and 
Coron, N, (1975) Astrophys . 210, L45, 

Barlow, M.O., and Silk, 0. (1976) Astrophys . <J. 207, 131. 

* 

Berg! and. G.D, (1969) IEEE Spectrum 6, 41. 

Blitz, L., Fich, M., and Kulkarni, S. (1983) Science 220. 1233. 

Boyd. R.W. (1982) Infrared Phys . 22, 157. 


196 


Bratt, P.R, (1977) In Semf conductors and Semimetals . R.K. Willardson and 
A.C. Beer, Eds. (Hew Yorks Academic Press), Vol. 12, p. 39. 

Burton, W.B., and Gordon, M.A. (1978) Annu . Rev. Astron . Astrophys , 14, 
275, 

Burton, W.B., Gordon, M.A., Banla, T.M., and Lockman, F.J. (1975) 
Astrophys . 203, 30. 

Cooley, J.W., and Tukey, d.W. (1965) Hath , Comput. 19, 297. 

Cravens, T.E., and Dalgarno, A, (1978) Astrophys . 0,, 219, 750, 

Crawford, M.K., Genzel, R,, Townes, C.H., and Watson, D.H, (1985) 
Submitted to Astrophys . J, ■ 

Dalgarno, A,, and McCray, R. (1972) Annu . Rev . Astron . Astrophys . 10, 
375. 

deJong, T», Dalgarno, A., and Boland, W, (1980) Astron . Astrophys . 91, 

68 , 

Downes, D., Wilson, T.L., Beiging, J., and Wink, J. (1980) Astron . 
Astrophys . Sup pi . Ser . 40, 379. 

Draine, B.T. ( 197 8 ) Ast rophys . J. Supp 1 . S er . 36, 595. 


Encrenaz, Th,, Gautier, D,, Vapillon, U, and Verdet, J.P. (1971) 

Astron, Astrophys , 11, 431, 

* 

Fellgett, P,B, (1951) Doctoral Thesis, Cambridge University, Cambridge, 
England, 

Flower, D.R., Launay, J.M,, and Roueff, E, (1977) In Les Spectres des 
Molecules Simples au Laboratoire et an Astrophysigue (Li&ge: University 
de Li age Institut a'Astrophysique), p, 137, 

Forrest, W,J,, McCarthy, O.F,, and Houck, J.R* (1980) Astrophys , tJ, 240, 
L37, ... 

Furniss, I., Jennings, R.E., King, K.J., Lightfoot, 0,F., Emery, R*J», 
Fitton, B,, and Naylo, DiA* (1983) Mon, Not , R, Astron , Soc , 202, 859. 

Gatley, I», Becklin, E.E., Sellgren, K», and Werner, M,W. (1979) 
Astrophys . J, 233, 575, % 

Genzel, R,, Reid, M.O., Moran, J,M., Downes, D., and Ho, P.T.P. (1982) 

In “Symposium on the Orion Nebula to Honor Henry Draper", Ann , N_,)\ 

Acad . Sci . 395, 142, 

Gerola, H*, and Glass gold, A.E. (1978) Astrophys , J, Suppl , Ser. 37, 1. 
Gillett, F.C., Dereniak, E.L., and Joyce, R.R. (1977) Opt . Eng . 16, 544. 


Gillett, F,C,, Forrest, W.J., and Merrill, K*M» (1973) Astrophys, d. 


183, 87, 


Glspert, R,, Puget, 0,L,, and Serra, G, (1982) Astro n. Astrophys , 106, 
293, 

Glassgold, A,E*, Huggins, P.d,, and Schucking, E,L., Eds, (1982) 
"Symposium on the Orion Nebula to Honor Henry Draper", Arm, N,‘X» Acad , 
Sci , 395- 

Glassgold, A,E, and Langer, W,D, (1973) Astrophys , J.» 186, 859, 
Goldsmith, P.F, and Langer, W,D, (1978) Astrophys , J. 222, 881, 

Hall, d,d. (1962) Phys, Rev, 128, 68, 

Haller, E,E», Hueschen, M,R,, and Richards, P.L, (1979) Appl, Phys . 
Lett, 34, 495, 

Harding, D,S,, and Harding, A, K, (1982) Astrophys , J.* 257, 603, 

Hart, L, and Pedlar, A, (1976) Mon , Not , _R‘ Astron , Soc, 176, 547, 

Harwit, M> (1984) In Galactic and Extragal actic Infrared Spectroscopy 
(Dordrecht: D, Rei del ) , p. 145. 


200 

Harwit, M., Kurtz, N.T., Russell, R.W., and Smyers, $• (1981) Appl . 
Opt , 20, 3792, 


Harwit, H# and Sloan, N.O.A. (1979) Hadamard Transform Optics (New York: 
Academic Press), 

Hayes, H.A. and Nussbaumer, H. (1984) Astron . Astrophys . 134, 193, 

Henderson, A.P., Jackson, P.D,, and Kerr, F.O, (1982) Astrophys , 263, 
116, 

Herter, T,, Briotta, D.A. Or., Gull, G.E., Shure, H.A,, and Houck, J.R, 
(1982) Astrophys . J,. 259, L1Q9, 

Houck, J.R. and Briotta, D»A. Or. (1982) Infrared Phys , 22, 215, 

Houck, J.R, and Ward, D.B, (1979) Pub! . Astron . Soc , Pac , 91, 140, 

Isobe, S, (1980) PubT , Astron . Soc . Jpn . 32, 423. 

Jaffa, D.T., Harris, A. I., Silber, M., Genzel, R», and Betz, A.L. (1985) 
"Detection of the 370 urn 3 p 2 ~ 3 ?i Fine-Structure Line of Cl," Submitted 
to Astrophys , 0_. ( Letters ). 

Oaffe, D.T. and Pankonin, V, (1978) Astrophys , J, 226, 869. 

Jenkins, E.B, (1978) Astrophys . J. 220, 107. 


Jenkins, G.M., and Watts, D.G, (1968) Spectral Analysis and its 
Applications (San Francisco; Holden-Day). 

Johnson, H,M. (I960) In Nebulae and Interstellar Hatter , B.M. 
Middlehurst and L*H. Aller, Eds. (Chicago; Univ. Chicago Press), p. 65, 

Kafatos, M*, and Lynch, J.P. (1980) Astrophys , J,. Suppl . Ser, 42, 611, 

Kazanskii, A, 6,, Richards, P,L,, and Haller, E,E, (1977) Appl. Phys, 
Lett, 31, 496, 

Keyes, R,J,, and Quist, T,M» (1970) In Semiconductors and Semimetals, 
Vol, 5 (New York; Academic Press), p. 321, 

Krugel, E,, Thum, C., Martin-Pintado, J,, and Pankonin, V, (1982) 
Astron . Astrophys , Suppl , Ser , 48, 345, 

Kulkarni, S.R., Blitz, L, } and Halles, C, (1982) Astrophys . J, 259, L63 

Kurtz, N.T,, Smyers, S,D,, Russell, R,W», Harwit, M,, and Mel nick, G, 
(1983) Astrophys , J,. 264, 538, 

Lang, K,R, (1980) Astrophysical Formulae , 2nd Ed. (New York; Springer- 
Verlag) , : 

Langer, W, (1976) Astrophys , _d, 206, 699, 


202 


Laques, P. and Vidal, d.-L. (1979) Astron , Astrophys. 73* 97. 

Launay, d.M. and Roueff, E. (1977) Phys . B 10, 879. 

Leung, C.M. (1975) Astrophys. d. 199, 340. 

McCarthy, D.E. (1967) J,. Opt. Soc, Am. 57, 699. 

McKee, C.F. and Ostriker, d.P. (1977) Astrophys . J., 218, 148, 

Melnick, G., Gull, 6.E., Harwit, M,, and Ward, D.B, (1978) Astrophys, 
d. 222, L137, 

Melnick, G,, Gull, G.E., and Harwit, M» (1979) Astrophys , d,. 227, L29. 

Melnick, G,d. (1981) Ph.D, Dissertation, Cornell University, Ithaca, NV. 

Melnick, G.d., Stacey, G.d>, Viscuso, P*d», and Fuller, C.F, (1985) 
"Observations of the 157.7 m CCHO Emission from the Galactic HII 
Regions W3 and W51," Submitted to Astrophys. J_. 

Mendel! , W.W, and Low, F.d. (1970) J, Geophys , Res , 75, 3319, 

Mendoza, C. (1983) In Planetary Nebulae , D.R. Flower, Ed, (IAU Symp. 
#103), Ch. 4, p, 143, 


Mezger, P.G. (1978) Astron , Astrophys , 70, 555. 

Mezger, P.G, (1984) In Galactic and Extra galactic Infrared Spectroscopy 
(Dordrecht: D. Reidel), p, 423, 

Mirabel), I.F. and Morras, fU (1984) Astrophys . 279, 86, 

* 

Moller, K.D., and Rothschild, W.G. (1971) Far-Infrared Spectroscopy (New 
York: Wi ley-Interscience). 

Moorwood, A.F.M., Baluteau, d,-P„ Anderegg, M», Coron, N,, Biraud, Y., 
and Fitton, B» (1980a) Astrophys , 238, 565, 

Moorwood, A,F,M, } Salinari, P,, Furniss, I., Jennings, R.E., and King, 
K,J, (1980b) Astron, Astrophys . 90, 304, 


Morton, D,C, (1975) Astrophys , 197, 85, 

Morton, D.C, (1978) Astrophys , 222, 863, 

Nishimura, T., low, F.O,, and Kurtz, R,F, (1980) Astrophys , 0., 239, 

L101, 

Osterbrock, p.E. (1974) Astrophysics of Gaseous Nebulae (San Francisco: 
Freeman), 

Papoulis, A, (1977) Si gnal Analysis (New York: McGraw-Hil 1 ) . 


Pelmbert, M, (X982) In "Symposium on the Orion Nebula to Honor Henry 
Draper", Arm. N.Y, Acad. Scl. 395 s 24, 

Phillips, T.G., Huggins, P.d., Kulper, and Miller*' R.E. (1980) 

Astrophys . d, 238, L103. 

Pottasch, S.R. (1965) Vistas Astron . 6, 149. 

Pottasch, S.R., Belntema, D.A., Ralmond, E», Baud, B., van Duinen, R», 
Hating, H#d,, Houck, d.R., dedong, T,, dennings, R.E., Olnon, F.H., and 
Wessellus, P.R. (1984) Astronhys . d. 278, L33, 

Pottasch, S.R., Wessellus, P.R., and van Duinen, R.d, (1979) Astron . 
Astronhys, 74, L15. 

Price, P.d, (1961) Phys . Rev. 124, 713, 

Rybickl, G.B. and Llghtman, A.P. (1979) Rad1at1ve i Processes In 
Astrophysics (New York: Wiley and Sons), 

Russell, R.W., Mel nick , 13., Gull, G.E., and Harwlt, M. (1980) 

Astrophys. £» 240, L99. 

Russell, R.W., Mel nick, G,, Snyers, S.D., Kurtz, N.T., Gosnell, T.R., 
Harwlt, M„ and Werner, M.W. (1981) Astrophys . £.250, L35. 


205 


Salpeter, E.E. (1979) In The Large Seale Character* sties of the Galaxy 
(IAU Symp, m)> p. 245. 

Sanders, D#B, (1981) Doctoral Thesis, State Univ. of New York (Stony 
Brook), p. 140# 

Sanders, D.B., Solomon, P.H., and Scovllle, N,Z, (1984) Astrophys* 

275, 182. 


Sanders, P.B., Clemens, D.P., Scovllle, N»Z„ and Solomon, P.H. (1984a) 
"The Massachusetts-Stony Brook CO Survey", to appear in Astrophys . 
Suppl. Ser , 

Saraph, H.E. (1973) 5. Phys , B 6, L243, 

Schmid-Burgk, J. (1982) In Numerical Data and Functional Relationships 
( New Series ). Group VI, Vol. 2c, H.H. Undolt and R. Bornstein, Eds. 
(Springer-Verlag: Berlin), p. 116, 

Schmidt, M. (1957) Bull . Astroft, Inst, Neth , 13, 247. 

Schraml, J, and Hezger, P.G. (1969) Astrophys . \L> 156, 269, 

Scoville, N.Z. and Solomon, P.M. (1975) Astrophys, J, 199, 1105, 

Shure, Houck, J.R., Gull, G.E., and Herter, T, (1984) Astrophys. 
■0, 281, L29. . . ■ 


Shure, M.A,, Houck, d.R., Herfcer, T# f Gull, G#E„ and Graf, P. (1984a) 
"Mid-Infrared Spectroscopy of Planetary Nebula* 1 , In proceedings of; 
"NASA/AtS.P. Symposium on Airborne Astronony'*, NASA/Ames Research 
Center, Moffett Field, CA, 

* 

Silk, dt (1973) Ann , Rev, Astron, Astrophys , 11, 269, 

Simpson, d,A. (1975) Astron , Astrophys , 39, 43, 

Solomon, P,M» and Sanders, D,B, (1980) in Giant Molecular Clouds in the 
Galaxy^ P.M, Solomon and M,G, Edmunds, Eds, (New York: Pergamon), p, 65, 

Spitzer, L., dr, (1978) Physical Processe s in the Interstellar Medium 
(New York: Wiley and Sons), 

Spitzer, L,, dr, and Tomasko, M.G. (1968) Astrophys , jJ, 152, 971, 

Stacey, G,d,, Kurtz, N»T», Smyers, S,D,, Harwit, M», Russell, R,W», and 
Mel nick, G, (1982) Astrophys. d, 257, L37, 

Stacey, G,d., Kurtz, N.T», Snyers, S.D., and Harwit, M. (1983) Mon , 

Not , R:, Astron , Soc , 202, 25P, 

Stacey, G.d„ Snyers, S,D,, Kurtz, N.T., and Harwit, M* (1983a) 
Astrophys . d. 265, L7, 


207 


Stacey, G.J., Snyers, S.D., Kurt?, N.T., and Harwit, H* (1983b) 
Astrophys , 268, L99, 

Stacey, 8,0. , Viscuso, P.J., Fuller, C.E., and Harwit, M. (1985) 
Astrophys. 0, 289 (in press). 

Stephens, T.L, and Dalgarno, A. (1973) Astrophys . 0. 186, 165. 

Storey, J.W.V., Watson, D.M., and Townes, C.H. (1979) Astrophys . J. 
233, 109, 

Strong, A.W., Riley, P.A., Osborne, J,U, and Hurray, O.D* (1982) 
Mon . Hot . R. Astron . Soc. 201, 495# 

Helens* A.G.G.M,, and Hollenbach, D. (1985) "Photodissociation 
Regions; , I. Basic Model". Submitted to Astrophys . 0,. 

Helens, A.G.G.M., and Hollenbach, D. (1985a) "Photodissociation 
Regions; II. A Model for the Orion Photodissociation Region". 
Submitted to Astrophys, _J, 

Todd, T.R., Clayton, C.M., Telfair, W.B., McCubbin, T,K», Jr., and 
Pliva, d. (1976) J, Molec . Spectrosc . 63, 201, 

Traub, W.A. and Stier, M.T. (1976) Appl . Opt , 15, 364,. 


208 


Tucker, K.D., Kutner, H.U, and Thaddeus, P. (1973) Astrophys . 

180, L13, 

Vidal, J.«L. (1982) In "Symposium on the Orion Nebula to Honor Henry 
Draper", Ann. N.Y, Acad. Scl. 395, 176. 

Viscuso, P.J., Stacey, G.O., Fuller, C.E., Kurtz, N.T., and Harwlt, M» 
(1985) "Submillimeter Observations of OH and CH in Orion", submitted to 
Astrophys. J. 

Walmsley, C.M, (1975) In HII Regions and Related Topics, T.L, Wilson and 
D. Downes, Eds. (Nev/ York: Springer-Verlag), p. 17, 

Ward, D.B., Dennison, B,, Gull, G., and Harwit, M. (1975) Astrophys . 

J. 202, 131. 

Watson, D.H, (1982) Ph.D, Dissertation, Univ. California, Berkeley, CA. 

Watson, D.M., and Storey, d.W.V, (1980) Int. J,. Infrared Millimeter 
Waves 1, 609. 

Watson, W.D. (1972) Astrophys . J.. 176, 103. 

Werner, M.W. (1970) Astrophys . Lett . 6, 81. 

Werner, M.W. (1982) In "Symposium on the Orion Nebula to Honor Henry 
Draper", Ann. N.Y. Acad. Sci . 395, 79, 


209 


Werner, M.W., Dlnnersteln, H.L., Hollenback, D.O,, Lester, D.F., Genzel, 
R., Watson, D.M., and Storey, J.W.V. (X982) Bull . Am, Astron. Soc. 14, 
610, 

Werner, M.W,, Gatley, I#, Harper, D.A., Beck! in, E.E., Lowenstein, R.F., 
Telesco, C.M., and Thronson, H,A, (1976) Astrophys , 204, 420. 

Westerhout, G, (1957) Bull . Astron. Inst . Neth . 13, 201, 

Westerhout, G, (1958) Bull . Astron . Inst. Neth. 14, 215. 

Wright, E.L. (1976) Astrophys , J,. 210, 250, 

Wynn-Wmiams, C.G., Genzel, R., Becklln, E.E., and Downes, D» (1984) 
Astrophys. J. 281, 172, 

York, D.G., Spltzer, L,, Bohlin, R.C., Hill, J,, Jenkins, E.B., Savage, 
B.D., and Snow, T.P, (1983) Astrophys . _J» 266, L55, 


